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DIFFUSION APPROXIMATION FOR A PROCESSOR SHARING 
QUEUE IN HEAVY TRAFFIC 1 

By H. Christian Gromoll 
EURANDOM 

Consider a single server queue with renewal arrivals and i.i.d. 
service times in which the server operates under a processor shar- 
ing service discipline. To describe the evolution of this system, we 
use a measure valued process that keeps track of the residual ser- 
vice times of all jobs in the system at any given time. From this 
measure valued process, one can recover the traditional performance 
processes, including queue length and workload. We show that under 
mild assumptions, including standard heavy traffic assumptions, the 
(suitably rescaled) measure valued processes corresponding to a se- 
quence of processor sharing queues converge in distribution to a mea- 
sure valued diffusion process. The limiting process is characterized as 
the image under an appropriate lifting map, of a one-dimensional re- 
flected Brownian motion. As an immediate consequence, one obtains 
a diffusion approximation for the queue length process of a processor 
sharing queue. 

1. Introduction. Consider a queueing system which consists of a single 
server with an infinite capacity buffer, to which jobs arrive according to a de- 
layed renewal process. The ith such arrival requires an amount of processing 
time that is the ith member of a sequence of independent and identically 
distributed strictly positive random variables. The server, rather than pro- 
viding service to just one job at a time, operates under a processor sharing 
discipline; that is, it works simultaneously on all jobs currently in the buffer, 
providing an equal fraction of its attention to each. Thus, at any given time 
that the buffer is nonempty, each job in the buffer is being processed at 
a rate that is the reciprocal of the number of jobs in the buffer. When the 
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server has fulfilled a given job's service time requirement, the job exits the 
buffer. This system is known as a processor sharing queue. 

The processor sharing service discipline can be viewed as an idealization 
of a round-robin or time-sharing protocol used in computer and communica- 
tion systems. There is a considerable literature on processor sharing queues 
(see [18] for a survey up to 1987), much of which assumes either Poisson 
arrivals or exponential service times. For a discussion of more recent work, 
including a handful of results for the case of generally distributed interarrival 
and service times (the GI/GI/l processor sharing queue), see [6]. 

In this paper, we present a heavy traffic diffusion approximation for a mea- 
sure valued process that keeps track of the "state" of a GI/GI/l processor 
sharing queue. A direct consequence of this is the existence of a diffusion 
approximation for the queue length process. Note that since the workload 
process in a GI/GI/l queue is the same for all nonidling service disciplines, 
the heavy traffic approximation for the workload process under a proces- 
sor sharing service discipline is the same as the well-known approximation 
under a FIFO (first-in-first-out) service discipline [9]. However, this simple 
relationship does not hold for the queue length process. 

The measure valued process that we study keeps track of the residual ser- 
vice times of jobs in the buffer. The residual service time at time t > of a job 
which has entered the buffer by time t, is given by the amount of processing 
time originally requested by the job minus the total amount of processing 
time it has received by time t. Jobs with residual service times at time t equal 
to zero have received enough processing time to fulfill their requirement and 
have departed the buffer. Let Mf denote the space of finite, nonnegative 
Borel measures on M + = [0, oo). The measure valued process : t > 0} is 
such that for each t > 0, fi(t) is the random element of A4p that has a unit of 
mass at the residual service time of each job currently in the buffer at time 
t. From this process, one can recover information about the performance of 
the system. For example, let Z(t) denote the number of jobs in the buffer, 
or queue length, at time t > 0. Then Z(t) can be recovered as the integral 
of fi{t) against the function that is identically one: Z(t) = {1, [/,(£)), where 
(1, //(£)) = / R+ lfj,(t)(dx). Similarly, let W(t) denote the sum of the resid- 
ual service times of all jobs in the buffer, or workload, at time t > 0. Then 
W(t) = (x,fi(t)), where = /k+ x{v)n{t)(dx) and x(x) = x. Although 

the process /u(-) includes information about the queue length and workload 
processes, it also provides a more detailed description of the state of the 
system than is available from these one-dimensional performance processes 
alone. It is this level of detail which facilitates our analysis. The process /x(-) 
is called the state descriptor for the processor sharing queue. Note that this 
terminology is not intended to imply that //(•) is necessarily a Markovian 
state descriptor (it may not be, since we do not include the residual interar- 
rival time, the time remaining until the next job arrives to the buffer, in the 
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state description). The process //(•) has previously been used by Grishechkin 
[5], along with other measure valued descriptors, in his heavy traffic analy- 
sis of the steady state distribution of a processor sharing queue. It has also 
been used recently by Gromoll, Puha and Williams [6], Puha and Williams 
[16] and Puha, Stolyar and Williams [15] for obtaining fluid limit results. 
Similar measure valued descriptors have also been used recently to describe 
other queueing systems. Doytchinov, Lehoczky and Shreve [3] and Kruk, 
Lehoczky, Shreve and Yeung [11] have used measure valued processes in 
the context of a queueing system with deadlines. Limic [12, 13] has used 
measure valued processes in studying the heavy traffic behavior of a LIFO 
(last-in-first-out) preemptive resume queue. 

In [6], a fluid (or law of large numbers) approximation for the state de- 
scriptor n(-) of a heavily loaded processor sharing queue was studied. It was 
conjectured there that an understanding of the steady state behavior of this 
fluid approximation would provide a key ingredient for establishing a heavy 
traffic diffusion approximation for the process //(■)• This is verified in the 
present paper and the role played by fluid approximations is informally de- 
scribed in the following paragraphs. The reader is referred to Sections 3 and 
4 for the complete treatment. 

For each r > in some sequence which tends to infinity, the fluid scaled 
and diffusion scaled versions of /i(-), denoted fi r {-) and jl r (-), are defined 
by pl {t) = r _1 /i(rt) and ff{t) = r _1 / u(r 2 t) for t > 0. Since (f{t) = j2 r (rt) for 
each t > 0, one can think of the diffusion scaled process p, r (-) over a finite 
time interval [0,T], where T > 1, as corresponding to the fluid scaled process 
fL r (-) over the time interval [0,rT]. The latter interval grows without bound 
as r — > co. To study Jl r {-) over this "long" time interval, one considers it in 
sections by covering [0, rT] with the overlapping finite time intervals [m, m + 
L], where m = 0, . . . , \ rT\ and L > 1 is fixed (here |_'J denotes the integer 
part). The shifted fluid scaled processes p, r ' m (-) are then defined for each m < 
\rT\ and t £ [0, L] by p, r ' m (t) = fx r {m + t). Thus, we can study the diffusion 
scaled process /F(-) over [0,T] by studying the family of shifted fluid scaled 
processes {fl r ' m (•) : m < [rT\} over [0, L\. By building on the techniques 
developed in [6], it is shown that as r^co, "good" sample paths of the 
processes {fl r ' m (-),m < \rT\} can be uniformly approximated on [0,L] by 
measure valued functions ((■) : [0, co) — ► My known as fluid model solutions 
(see Definition 4.2 in Section 4). By "good" sample paths, we mean sample 
paths in a set whose probability approaches 1 as r — ► co. 

Recently, Puha and Williams [16] have shown that under mild conditions, 
such fluid model solutions £(•) converge to a steady state; that is, for any such 
((■), there is a Coo G -M.F such that as t — > co, ((t) — ► (^0 in the topology of 
weak convergence of measures. Moreover, Coo depends only on the limiting 
service time distribution v, and on the first moment of C(t)i which remains 
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constant in t. Indeed, 

Coo = A„(x,C(<)) for alii >0, 

where A u : R + — ► A^f is the lifting map defined by A u w = w{x, v e )~ l v e and 
v e denotes the residual excess lifetime distribution associated with is; that is, 
u e is the probability measure on R + satisfying v e ([0, x]) = (%, v) -1 u{y, oo) dy 
for all x G R+. In [16], conditions are also identified under which convergence 
to steady state occurs uniformly for certain collections of fluid model solu- 
tions. 

In the present paper, the steady state results in [16] are combined with 
the fluid approximations of {p, r ' m (■) : m < [rT\} described above to show 
that if L > 1 is sufficiently large and t 6 [L — 1, L] , then for "good" sample 
paths, 

(1.1) £ r,m (t)«A„<x,£ P ' m (t)). 

This fact is applied to the diffusion scaled state descriptors in the following 
way. For any t € [(L — l)r ,T], there is an m £ {0, . . . , L r ^J} an d an s € 
[L — 1,L] such that 

(1.2) (i r {t) = fl r (rt) = fL r (m + s) = U r ' m {s). 

Combining (1.2) with the approximation (1.1) yields that for "good" sample 
paths, and for t € [(L — l)r _1 ,T], 

(1.3) A r (t)«A„< X ,£'(t)). 

It is shown that under mild conditions, including standard heavy traffic 
assumptions, the above approximation is valid on the entire interval [0,T], 
with probability approaching 1 as r — > oo. Recall that (x, /* r (^)) = r _1 W (r 2 t) 
is the (diffusion scaled) workload at time t, which we denote by W r (t). Thus, 
the essence of (1.3) is that as the diffusion scaled process /i r (-) approaches 
the heavy traffic limit, it can be recovered from the diffusion scaled workload 
process W r (-) on [0, T] by an appropriate lifting map. This phenomenon is 
known as state space collapse. Since it is known that under the assumptions 
we will impose, the process W r (-) converges in distribution to a reflected 
Brownian motion W*(-) on K + , a diffusion approximation for (i r {-) follows 
quickly from state space collapse. (See Sections 3 and 4 for more details.) 
This yields a limiting measure valued diffusion process which is confined to 
the one-dimensional subspace {cv e : c > 0} of A^p. It is given by 

fi*(-) = A u W*(-). 

A consequence of the above diffusion approximation for the diffusion 
scaled state descriptor fi r (-) is the existence of a diffusion approximation 
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for the diffusion scaled queue length process Z r (-) = (l,/x r (-)). The limiting 
process is the one-dimensional reflected Brownian motion 

Z*(.) = <l,//(-)>. 

The precise form of Z*(-) is given in Corollary 2.4, which verifies a conjecture 
of Harrison and Williams [8] in the case of a single processor sharing queue 
with a single class of jobs. We note that in cases where the process Z*{-) has 
negative drift, the steady state distribution of Z*(t) as t — > oo is consistent 
with that obtained by Grishechkin [5] using other means (see Section 2.4). 

The method outlined above for obtaining a heavy traffic diffusion approx- 
imation using fluid approximations as an intermediate step is analogous to 
the method developed by Bramson [2] and Williams [17] for proving diffu- 
sion approximations for open multiclass networks with HL (head-of-the-line) 
service disciplines. Indeed, the framework developed in [2, 17] serves as the 
primary motivation for the approach taken here for the processor sharing 
discipline. Note that since processor sharing is not an HL service discipline, 
a large part of the machinery needs to be developed from first principles. 

The paper is organized as follows. In Section 2 we give a precise description 
of our model and assumptions, and we state our main results. Section 3 is 
devoted to fluid scale analysis, which prepares us for the proof of state space 
collapse. The latter is given in Section 4. 

1.1. Notation. The following notation will be used throughout the paper. 
Let N = {1, 2, . . . } and let R denote the set of real numbers. For a, b € M., we 
write a V b for the maximum of a and b, a A b for the minimum of a and b, a + 
for the positive part of a, [a\ for the largest integer less than or equal to a, 
and \a\ for the smallest integer greater than or equal to a. The nonnegative 
real numbers [0, oo) will be denoted by R+. For a function <?:M+ — >M, let 
llfflloo = sup a . GR+ \g(x)\ and \\g\\ K = sup^^ \g(x)\ for each K > 0. 

For a set B C R+ , we denote the indicator of the set B by 1b • We also 
define the following real valued functions on R + : x( x ) = x i f° r x € R+, and 
<p(x) = 1/x, for x G (0, oo) with ip(0) = 0. For a topological space A, denote 
by Cb(A) the set of continuous, bounded, real- valued functions defined on 
A. In addition, for an interval JcR, C^(J) is the set of once continuously 
differentiable, real-valued functions defined on / that together with their 
first derivatives are bounded on I. 

Recall that Aip is the set of finite, nonnegative Borel measures on R + . 
Consider £ € A^f an d a Borel measurable function — ► R which is 
integrable with respect to £. We define {g, (} = J R+ g(x)((dx). Our equations 
will involve expressions of the form Jj a ^ g(x — a)£(dx) for a > 0. To ease 
notation throughout, we write this as (g(- — a),Q, making the convention 
that such a g is always extended to be identically zero on (-co, 0). The space 
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Mf is endowed with the topology of weak convergence of measures; that 
is, for Cn.,C € Mp, n G N, we have Q n — > ( if and only if (g,(n) — ► {QiO 
as n — > oo for all 5 : R_|_ — ► R that are bounded and continuous. With this 
topology, M.f is a Polish space [14]. It will be convenient to deal with weak 
convergence on Mp by means of a suitable metric which we now define. First, 
choose a countable set of nonnegative functions {g G }f° =1 C C£(R+) that are 
each bounded by 1, and such that {g G is convergence determining for 
A4f, that is, such that Q n -^-> ( in M-f, as n — > 00, if and only if for each 
fceN, (fff , Cn) — (fff 5 0) as n ^ 00. We refer the reader to [4], Chapter 3, 
Section 4, for a possible construction. Note that the construction in the proof 
of Proposition 4.2 in [4], Chapter 3, Section 4, can be modified to require that 
g G k G C^(R+) and ||fl^||oo < 1 for k G N. Next, for each k G N, let h G k G Cj(M+) 
be a nonnegative function that is bounded by 1 and satisfies Jvl = on 
[0, fcg = 1 on [fc.oo) and ||(^)'||oo < 2. Letting ^ = {^}^ =1 U{^}^ =1 , 

we define a metric on A4f in terms of Q as follows. For Ci 5 C2 £ -^f, define 

00 

d[Ci,0J = E2-*(|{rf > Ci>-(rf,0»>|Ai) 

(L4) +sup|(^,Ci)-(^,C 2 )|. 

fceN 

It is straightforward to see that d[-,-] is a complete metric on A4p that 
induces the topology of weak convergence of measures. Note that although 
d[-, •] depends on the set Q, we do not indicate this explicitly in the notation. 
We denote the zero measure in M.f by and the measure in Mf that puts 
one unit of mass at the point x G R+ (or Dirac measure at x) by 5 X . We will 
use <5+ to denote the following truncated Dirac measure: 



r+ f&z, X G (0,00), 

x \0, x = 0. 



We use "^=^" to denote convergence in distribution of random elements of 
a metric space. For random elements X, Y of a metric space, X ~ Y will de- 
note equivalence in distribution. Following Billingsley [1], we will use P and 
E, respectively, to denote the probability measure and expectation operator 
associated with whatever space the relevant random element is defined on. 
All stochastic processes used in this paper will be assumed to have paths that 
are right continuous with finite left limits (r.c.1.1.). For L G (1, 00) and a Pol- 
ish space 5, we denote by T)l(S) = D([0, L],S) [resp. Doo(«S) = D([0, 00), S)] 
the space of r.c.1.1. functions from [0,L] (resp. [0,oo)) into 5, and we endow 

this space with the usual Skorohod Ji-topology [4]. We use "-—!->" to denote 
convergence in this topology. 
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2. The processor sharing queue. In this section we give a precise descrip- 
tion of the processor sharing queue, specify our assumptions and state our 
main result. A formal definition of the processor sharing queue, as considered 
in this paper, was previously given in [6]. The reader is referred there for 
a detailed discussion. However, since the diffusion approximation presented 
here requires stronger assumptions on the model than were made in [6] , it is 
necessary to briefly review the definition. This is the subject of Section 2.1, 
where we introduce a sequence of processor sharing queueing models and the 
associated notation. Section 2.2 describes the scaling and time shifts we will 
apply to this sequence of models, and Section 2.3 specifies our asymptotic 
assumptions on the sequence. The statement of our main result appears in 
Section 2.4. 

2.1. A sequence of processor sharing queues. We now specify a sequence 
of processor sharing queueing models indexed by r G 1Z, where 1Z C (0, oo) 
is a sequence which increases to infinity. Each model in the sequence may 
be defined on a separate probability space and we use P r and E r , respec- 
tively, for the probability and expectation operator on each of these spaces. 
The rth model in the sequence consists of the following: A server, which 
processes jobs from an infinite capacity buffer according to the processor 
sharing discipline, a pair of stochastic primitive processes E r (-),V. T , which 
describe the arrival of work to the buffer, and a random initial condition, 
which specifies the state of the system at time 0. A measure valued process 
//"(") together with a set of descriptive equations describe the time evolution 
of the state of the system. 

The exogenous arrival process E r {-) is a rate a r delayed renewal process 
associated with a sequence {u'^}f^ 1 of finite nonnegative interarrival times. 
For t > 0, 

E r (t) represents the total number of jobs which have arrived to the buffer 
during the time interval (0, t] . We refer to the ith job to arrive to the buffer 
after time simply as the ith job. In contrast, jobs which are already in the 
buffer at time will always be referred to as initial jobs. The quantity u\ is 
the arrival time of the first job, and u^, i > 2, is the elapsed time between 
the arrival of the (i — l)st and ith jobs. Thus, for i > 1, U[ = J2)=i M j 1S t ne 
arrival time of the ith job. We define Uq = 0. So for t > 0, 

(2.1) E r (t)=svLp{i>0:Ul <t}. 

It is assumed that {u\}'^ l is a sequence of independent random variables and 
that {it[}^ 2 is i.i.d. with mean (a r )~ 1 G (0, oo) and standard deviation a r < 
oo. The first element u\ of the sequence is assumed to be strictly positive 
with finite mean. 
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The service process {Vf,i = 1,2,...} records the cumulative amount of 
processing time required from the server by the first i jobs. It is defined 
from a sequence of strictly positive service times by 

i 

(2.2) V[ = J2^ for i = l,2,.... 

3=1 

The quantity v\ represents the amount of processing time that the ith job 
requires from the server. It is assumed that {w i r }^ 1 is a sequence of strictly 
positive i.i.d. random variables with common distribution given by a Borcl 
probability measure v T on R + . We assume v r has mean (/? r ) _1 £ (0,oo) and 
standard deviation b r < oo. Finally, we assume that the service process V. r 
is independent of the exogenous arrival process E r (-). 

The initial condition specifies Z r (0), the number of initial jobs present 
in the buffer at time zero, as well as the service time requirement for each 
of these initial jobs. We assume Z r (0) is a nonnegative, integer valued ran- 
dom variable. The service times for initial jobs are taken to be the first 
Z r (0) elements of a sequence {v"j}'i£ =1 of strictly positive random variables. 
The random variables Z r (0) and {vj}'jL 1 are assumed to be independent of 
{n[}^ 2 an d {vDi^i, but are not assumed to be independent of one another. 
A convenient way to express the initial condition is to define an initial ran- 
dom measure fi r (0), which is the nonnegative, finite random Borel measure 
on R + given by 

z r (o) 

(2.3) /i r (0)=£*5- 

3=1 

Recall that 5 X denotes the Dirac measure at x € R+. Henceforth, ^ r (0) will 
be used as the initial condition. We assume that fi r (0) satisfies 

(2.4) El<l,//(0))]<oc, 

(2.5) El< X ,//(0)>]<oo. 

Note that since (l,//(0)) = Z r (0) and (x,//(0)) =Y%Li*Vj, assump- 
tions (2.4) and (2.5) mean, respectively, that the expected initial queue 
length and expected initial workload are finite. 

The state descriptor is a performance process which describes the time 
evolution of the state of the processor sharing queue. It is a measure valued 
process /i r (-) : [0, oo) — > M.f such that for each t > 0, the random Borel 
measure /i r (t) has one unit of mass located at the residual service time of 
each job that is in the buffer at time t. The residual service time at time 
t > of job % < E T {t) [reap., of initial job j < Z r (0)\, denoted JSJ"(t) [resp., 
Rj(t)], is the remaining amount of processing time required to fulfill the 
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service time requirement of the job. If this residual service time is zero, the 
job has completed service and has departed the buffer. Thus at t > 0, the 
state descriptor is given by 

Z r (0) E r (t) 

( 2 -6) /*"(*)= E E 

j=l 3 i=l 

Recall that 6% is the Dirac measure at x for x € (0, oo) with 5$ = 0, which 
ensures in the above that has one unit of mass only for each job with 

nonzero residual service time at time t > 0. Let Z r (t) denote the number 
of jobs in the buffer, or queue length, at time t > 0. Then clearly, Z r (t) = 
(1, //(£)) for all t > 0. Thus, under the processor sharing discipline, any job 
that is present in the buffer at time t receives service at the instantaneous 
rate (l,// r (t)) -1 . Note that if a job is present in the buffer at time t, then 
(1, //(£)} ^ 0. The cumulative service per job provided by the server up to 
time t > is defined by 

(2.7) S r (t)= [\{{l^ r {s)))ds, 

Jo 

where <p(x) = 1/x for x E (0, oo) with (p(0) = 0. Then, since job % < E r (t) 
arrives at time XJ\ < t, the cumulative amount of processing time that job i 
receives by time t is equal to v\ A (S r (t) — S r (U()). Similarly, the cumulative 
amount of processing time that an initial job j < Z r (0) receives by time t is 
equal to Vj A S r (t). Therefore, for t > 0, the residual service times are given 
by the equations 

(2.8) R\{t) = « - S r (t) + S r m)) + for i = 1, . . . , E r (t), 

(2.9) R^t) = ($ - S r (t)) + for j = 1, . . . , Z r (0). 

Given the primitive processes E r (-),V. r , and the initial condition /i r (0), 
equations (2.6)-(2.9) determine the state descriptor //(•), the process S r (-) 
and the residual service times. This fact is not difficult, although somewhat 
tedious, to show. We note that the definition of the state descriptor given 
here is equivalent to the formulation in [6]. 

Let W r (t) denote the (immediate) workload in the buffer at time t > 0. 
This is defined as the amount of time the server would have to work in order 
to complete the remaining service time requirements of all jobs in the buffer 
at time t, assuming no new arrivals take place. Since this is given by the 
sum of the residual service times of all jobs which are in the buffer at time 
t, we have 

(2.10) W r (t) = { X ,fi r (t)) fort>0. 

The cumulative service per job process S r (-) will play a particularly im- 
portant role in our analysis. We will find it convenient to have notation for 
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the increments of this process. For t,h>0, define the cumulative service per 
job in [t, t + h] by 

rt+h 

(2.11) Sl t+h = S r (t + h)-S r (t) = jf ^«l,//( S )»d S . 

Then for UJ <t, the amount of service received by the ith job by time t 
can be written as v\ A S T ur t , and the residual service time at time t can be 
written 

Rl(t) = (vI-S r UT ,t) + - 

2.2. Scaling. Our result concerns the asymptotic behavior of processor 
sharing queues on diffusion scale. In particular, the focus of our attention is 
the diffusion scaled state descriptor, which is defined for t > by 

(2.12) A r (*) = V(r a t). 

r 

We will also be interested in diffusion scaled versions of the workload and 
queue length processes, defined for t > by W r (t) = r~ 1 W r (r 2 t) = (x,fi r (t)} 
and Z r (t) = r Z r (r t) = (l,/2 r (t)), respectively. Much of our understanding 
of the diffusion scaled state descriptor /F(-) will be derived from results 
about the fluid scaled state descriptor, which is defined for t > by 

(2.13) J?(t) = V(rt). 

r 

The relationship fi r {t) = Jl r {rt) is essential for bootstrapping results from 
fluid scale up to diffusion scale. Recall that our approach for studying the 
diffusion scaled process £i r (-) over a fixed finite time interval [0, T], for T > 1, 
is to look at the fluid scaled process fi r (-) over the time interval [0,rT]. 
Specifically, we study overlapping sections of fl r (-), each defined on a finite 
time interval of fixed length L > 1. For each r S 1Z, t > and m G {0, 1, . . . }, 
define 

(2.14) fi r ' m (t)=fd r {m + t). 

Then since for each r G 1Z, the time interval [0, rT] is covered by the \rT\ + 1 
overlapping time intervals [m,m + L], for < m < \rT\, we have that for 
any t G [0, rT], there is (at least one) m < [rTj and s G [0, L] such that 

(2.15) fl r (t)=^ m (s). 

Much of Section 3 is devoted to extending results in [6] on the asymptotics 
(as r — > oo) of p, r (-) over [0, L], so that they hold for fi r ' m {-) over [0,L] for 
all m < [rT\ . For this discussion, we define the following fluid scaled and 
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shifted fluid scaled processes. For all r € 1Z, t > 0, h > and m G {0, 1, . . . }, 
let 



(2.16) 


E r {t) = 




E r ' m {t) 


= E r (m + t), 


(2.17) 


z r (t) = 


-Z r (rt), 
r 


Z r > m (t) 


= Z r (m + t), 


(2.18) 


W r (t) = 


-W r (rt), 
r 


W r ' m (t) 


= W r (m + t), 


(2.19) 


or 

°t,t+h — 


or 

D rt,r{t+h)i 


°t,t+h 


or 

— D m+t,m.+t+h 


Note that by (2.19), 









S t.?lh = S r(m+t),r( m +t+h) = I ' <P((1, (*))) ds 



\t+h u r{m+t),r(m+t+h) 

(2.20) 



r{m+t+h) 
r(m+t) 



t+h 

<p((i,pr> m (s)))d8. 



2.3. Heavy traffic assumptions. In this section, we specify the assump- 
tions under which the diffusion approximation will be proved. Let a > 0, 
a > and 6 > be fixed constants and let v be a probability measure on M + 
that does not charge the origin, and satisfies 

(2.21) < X 4 +V)<oc. 

Then j3 = (x, v) ~ 1 is positive and finite and b = ((x 2 , v) — (x, u) 2 ) 1 ! 2 is finite. 
In order to obtain convergence in distribution of the diffusion scaled state 
descriptors jl r {-) to a diffusion process, we impose the following asymptotic 
assumptions on the sequence of processor sharing queues defined in Sec- 
tion 2.1. For the sequence of arrival processes, we assume that as r — > oo, 

(2.22) (a r ,a r ) — >(a,a), 

(2.23) E r K]/r^0, 

(2.24) limsupE r [(n r 2 ) 2+e ] < oo. 

r— »oo 

For the sequence of service processes, we assume that as r — > oo, 

(2.25) v r ^ v, 

(2.26) (jr,b r )^ (P,b), 

(2.27) limsup(x 4+e ,^ r } <oo. 

r— *oo 

Define the traffic intensity parameter for the rth system by p r = a r / (3 r , and 
let p = a/f3. In addition to the above, our assumption that the sequence of 
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processor sharing queues approaches heavy traffic requires that 

(2.28) p = l. 

We also require the sequence to approach heavy traffic at an appropriate 
rate and assume that as r — > oo, 

(2.29) r(l - p r ) — ► A for some AeK. 

Assumption (2.23) implies that the initial residual interarrival time vanishes 
on diffusion scale. Assumptions (2.24) and (2.27) imply Lindeberg-type con- 
ditions, which along with (2.22), (2.23) and (2.26), are needed to imply func- 
tional central limit theorems for the triangular arrays {u\\i = 1, 2, . . . } T £ll 
and {yl\i = 1,2,... } r en anci ultimately, convergence in distribution of the 
diffusion scaled workload processes to a reflected Brownian motion (see 
Proposition 3.1). Assumption (2.27) is two moments stronger than what 
would normally be used for such functional central limit theorems. The ad- 
ditional restriction is used in a separate part of our analysis to estimate 
moments of the shifted fluid scaled state descriptors {/2 r ' m (-), m < [rT\} 
(see Lemma 3.6). 

We also make assumptions on the asymptotic behavior of the diffusion 
scaled initial measure £t r (0). Since fi r (0) = fi r (0), we make our assumptions 
in terms of £t r (0) since they will be used in that form. The following definition 
is central to these assumptions. 

Definition 2.1 (Invariant manifold). Let v e denote the excess lifetime 
distribution associated with u, that is, v e is the probability measure on Un- 
satisfying {l[ , x ],Ve) =PIo(^(y,oo)^) d y for all x£R + . Let 

Following usage in [2], we refer to the one parameter family of measures 
My C A4f as the invariant manifold associated with v. 

Let be a random measure taking values in Mp, such that 

(2.30) 6eM„ a.s., 

(2.31) E[(l,9)]<oo. 

Note that assumption (2.21) implies that {x,v e ) < oo and (x + , u e ) < oo, 
for k = 1,2,3. Combining this with (2.30) and (2.31) implies that also 
satisfies 



(2.32) 
(2.33) 



E[( X ,e)]<oo, 
E[(x fc+e ,©)] < oo for k = 1,2,3. 
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An important property of G that we will need is that it has no atoms a.s. 
This is a trivial consequence of (2.30) since v e has a density function. This 
"no atoms" property of will be used in the following form: 



(2.34) limP sup (1\ XX+K \,@) <e =1 for all e > 0. 

The equivalence of (2.34) to the fact that has no atoms a.s. is proved in 
[6], Lemma A.l. 

For the sequence of (scaled) initial measures /t r (0) = fl r (0), we assume 
that as r — > oo, 

(2.35) (fi r (o),(x,fi r (o)),(x l+e ,fi r (o))) (e,{ x ,Q),( x 1+e ,e)). 

Note that for any g € C 6 (R+), * fl : M F — > R defined by ^ g (Q = (g,() is 
a continuous function. So the first component of (2.35) implies by the contin- 
uous mapping theorem ([1], Theorem 5.1) that for any such g, (g,p, r (0)) 
(<?, 0) as r — > oo. The second component of (2.35) implies that the fluid/diffusion 
scaled initial workload converges in distribution, that is, V^ r (0) = W r (0) => 
(x, 0} as r — > oo. 

To simplify the statements of results for the remainder of the paper, we 
now summarize our assumptions in the following. 

(Q.l) There is a sequence of processor sharing queues defined in Section 2.1 
such that for some constants 6 > 0, a > 0, a > 0, some probability 
measure v on R + that does not charge the origin, and some random 
measure taking values in Mp, (2.21)-(2.31) and (2.35) hold. 

2.4. Main result. In order to state our main result, we will need the 
following definition. 

Definition 2.2. Assume (Q.l). Let A u :R + — ► Mf be the lifting map 
associated with v given by 

w 

(2.36) A v w = -. rv e forwGlR + . 

Theorem 2.3. Assume (Q.l). Then as r — > oo, the sequence of diffu- 
sion scaled state descriptors {fi r (-)} converges in distribution to the measure 
valued process = A U W*(-), where W*(-) is a reflected Brownian mo- 

tion on R + with drift —A, variance aa 2 + (5b 2 and initial value VK*(0) that 
is equal in distribution to (x>0)- 

Corollary 2.4. Assume (Q.l). Then as r — > oo, the sequence of dif- 
fusion scaled queue length processes {Z r (-)} converges in distribution to the 
process Z*{-) = C U W*(-), where W*(-) is the reflected Brownian motion spec- 
ified above, and C v = 2/3(1 + (i 2 b 2 )~ l . 
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PROOF. Since Z r (-) = (l,/t r (-)) and since — > R denned by 
^i(^) = is continuous, the result follows from Theorem 2.3 by the 

continuous mapping theorem. The form of the constant C v follows from 
(2.36) by rewriting the mean of v e in terms of the mean f}~ , and standard 
deviation b, of v. □ 

Harrison and Williams [8] have conjectured that for a queue operating 
under the processor sharing discipline (possibly with several different job 
classes which may arrive to the buffer) , a heavy traffic limit theorem should 
hold for the queue length process associated with each job class. They predict 
that for each job class, the limiting queue length process should be a constant 
multiple of the reflected Brownian motion obtained as the limiting workload 
process for the queue. Moreover, they specify the form of the constants for 
which the above claim should hold (see A. 58, A. 60 and A. 61 in [8]). Corollary 
2.4 verifies this conjecture in the case of a single job class. In particular, the 
constant C u validates (after some rewriting) the prediction made in [8], A. 60 
and A. 61. 

In cases where A is positive, the process Z*(-) has a steady state distri- 
bution which can be computed as in [7] to be 

lim P(Z*(t)<x) 

(2 37) 

= 1 - exp(-A(^- 2 + b 2 ){a 2 + b 2 Y l x). 

It is interesting to note that in the case /3 = 1, this is consistent with the 
steady state distribution obtained by Grishechkin [5]. Indeed, since the ap- 
proach taken in [5] was to first take the limit as t — > oo to get the steady 
state distribution of Z r {t) and then take the limit of the resulting steady 
state distributions as r — > oo, the consistency with (2.37) suggests that an 
interchange of limits is possible in this situation. 

The strategy for proving Theorem 2.3 can be divided into two main tasks. 
The first task is to establish a tightness property for the processes {/2 r,m (-), 
r G 7£, to < [rTj }. This is accomplished in Section 3. In Section 4, we under- 
take the second task, which is to identify limit points obtained in Section 3 
as fluid model solutions (see Definition 4.2) and to show that, consequently, 
the process ff{-) can be approximated, as r — > oo, by "overlapping" fluid 
model solutions. This will be done in such a way that as r — > oo, the random 
measure fi r (t) is close to the steady state of some fluid model solution for 
all t £ [0,T], with arbitrarily high probability. Combined with a result in 
[16] which provides rates of convergence to steady state for such fluid model 
solutions, this leads to the proof of state space collapse. Finally, the proof 
of Theorem 2.3 collecting these results appears at the end of Section 4. 
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3. Compactness. The aim of this section is to prove a tightness prop- 
erty for the shifted and fluid scaled state descriptors {fi r ' m (■) , r G 1Z, < m < 
[rT\ }. Specifically, our goal is to show in Corollary 3.16 at the end of this sec- 
tion that sequences indexed by r G 1Z of "good" sample paths of the processes 

r G 7£, to < [rT\} are relatively compact in ~Dl(A4f), for L > 1. Here, sam- 
ple paths are "good" with arbitrarily high probability as r — > oo. The road 
to Corollary 3.16 involves a detailed analysis of the behavior of the processor 
sharing queue on fluid scale and we rely heavily on the ideas developed in 
[6], Section 5. Nevertheless, many of the techniques developed there need to 
be refined and several new ideas are needed for the present analysis. This 
is to account for the fact that for each r G 1Z, we are dealing with order r 
measure valued processes (one for each m G {0, . . . , [^J}) instead of just 
one. 

The section is organized similarly to [6], Section 5. In Sections 3.1 and 3.2 
we describe a dynamic equation satisfied by the fluid scaled state descrip- 
tors {p> r ' m (-)} and a well-known heavy traffic limit theorem for the workload 
processes {(x>/" r,m ("))}> respectively. Section 3.3 contains a generalization of 
Lemma 5.3 in [6]. It states that under certain conditions, the fluid scaled 
queue length process (l,/x r (-)) can be bounded above over long time inter- 
vals. This will yield an upper bound for the shifted fluid scaled queue length 
process (l, / u r,m (-)) on [0,L] for each m < \rT\. Section 3.4 contains a func- 
tional weak law of large numbers estimate, which is a refinement of Lemma 
A.2 

in [6]. In Section 3.5, we obtain an upper bound for certain moments of 
the processes {yw r ' m (-)}- m Section 3.6, we state and prove Lemma 3.8. 
This lemma combines the results of Sections 3.1-3.5 for subsequent easy 
reference. Lastly, Section 3.7 contains further analysis of the fluid scaled 
processes {p- r ' m (-)}- It consists of four secondary lemmas which are conse- 
quences of the results summarized in Lemma 3.8. These four lemmas lead 
up to an oscillation bound (Theorem 3.14) for the processes {fi r ' m (-)} as 
well as Corollary 3.16. 

3.1. Dynamic equation. We begin by specifying a dynamic equation sat- 
isfied by the processes {A r,m (-)} f° r m — l r T\ - Starting with (2.6) and sub- 
stituting in the definition of the residual service times (2.8) and (2.9), one 
obtains, after some simplification, that for each r, a.s. for each Borel mea- 
surable function g : M + — ► R and all t, h > 0, 

(g,V r (t + h)) = ((t {0 ,oo)9)(--S r t>t+h )^ r (t)) 

( 3-1 ) E r {t+h) 

+ ( 1 (0,oo)5)« -Suv,t+h)- 

i=E r (t)+l 
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Recall that we always assume g is extended to be identically zero on (—00, 0) 
so that functions of the form g(- — a) are well defined on R + for any a > 0. 

Equation (3.1) takes the following form for the shifted fluid scaled pro- 
cesses £t r,m (-). For each r 6 TZ and each m < [rT\, a.s. for each Borel mea- 
surable function g : M + — ► R, and all t, h > 0, 

(g,fl r ' m (t + h)) = ((l ( o,oo)5)(- -S r ^ h W m {t)) 

v ' ; rE r < m (t+h) 

+ - E (%0,oc)S)K-^-i_ mit+/l ). 
i=r£ 7 '. m (t)+l 

We refer to (3.2) as the dynamic equation for j2 r ' m (-). Frequently we will set 
g = 1 in this equation, in which case it will look like 



l,fl r ' m (t + h)) = (t { sr :r+h ^fi r ' m (t)) 
1 



rE r < m (t+h) 



r 



■m,t+h' 



i=rE r > m {t)+l 



An important bound which we will use often is obtained from (3.3) by ig- 
noring any processing of jobs during the fluid scaled time interval [t,t + h]: 

(3.4) (1, i?> m {t + h)) < (l,ft r > m (t)) + E r > m (t + h)- E r > m (t). 

3.2. Diffusion limit for the workload process. We will take advantage of 
the following well-known result. 

Proposition 3.1. Assume (Q.l). Then as r^oo, the sequence of dif- 
fusion scaled workload processes \W r {-)} = {(x>A r ('))} converges in distri- 
bution to a process W*{-), which is a reflected Brownian motion on WL+ 
with drift —X, variance aa 2 + (5b 2 and initial value W*(0) that is equal in 
distribution to {x, @) ■ 

Proof. It is well known that assumptions (Q.l) are sufficient to imply 
the above result for the workload process of any single server, single class 
queue operating under a work conserving service discipline (including pro- 
cessor sharing) . The use of functional central limit theorems and continuous 
mappings to prove such results goes back to [9]. For a detailed account, 
see, for example, [17], which covers a much more general setting than that 
considered here. □ 



Corollary 3.2. Assume (Q.l) and let T > 1 and < rj' < 1 be given. 
Then: 
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(i) There exists M > 1 such that for any L > 1 , 

(3.5) liminfP r ( sup \\( X , ^ m (-))h < m) > 1 - V ' . 

(ii) For any L > 1 and any 7 > 0, 

(3.6) liminfW sup |<X,//' m (*)> - <X,//' m (0)>| < > 1 - rf. 

r ^°° \m<[rT\, te[0,L] 4 / 

Proof. By expanding the definition of /x r,m (-) and rewriting it in terms 
of diffusion scaling, we see that it suffices to show that: 

(i) There exists M > 1 such that for any L > 1, 



(3.7) liminfPM sup \W r (t)\ < M > 1 - 7/. 

r ~*°° \te{0,T+(L/r)] J 

(ii) For any L > 1 and any 7 > 0, 

(3.8) liminfP r ( sup \W r (t + h) - W r (t)\ < 1 ) > 1 - rf , 

r ^°° \te[o,T], he[0,L/r] 4 / 

both of which follow easily from Proposition 3.1, using the tightness of 
{W T (-)}r£ii and the fact that the limiting process W*(-) is a.s. continuous. 
Note that M can indeed be chosen independent of L. □ 

3.3. Stability of the queue length process. The following lemma is a sim- 
ple generalization of Lemma 5.3 in [6]. It will be used in two separate places 
to bound the total mass of the fluid scaled state descriptor over certain time 
intervals. 

Lemma 3.3. Assume (Q.l) and let T > 1 be given. Suppose that for 
a given r € 1Z, there is an event A r and there are constants to € [0,rT] and 
t\ , c, I > such that on A r , 

(i) SK Vt&[toM+tl] Er(t + l)-W(t)<l 

(ii) sup te[toA)+tl] {x,fl r (t)) < {, 

(iii) (l,/z r (to))<§. 

Then on A r , 



(3.9) 



sup (1, //(£)) < c. 
te[t ,to+*i] 
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Proof. Cover [to, to + 1\] with the time intervals Ik = [to + kl, to + (k + 
1)1], where k 6 {0, . . . , L^i/^J}- We prove by induction on k that on A r , 

(3.10) sup(l,£ r (*))<c. 

te/ fe 

We first verify the case k = 0. We have on ^4 r , for t £ Io, 
(l,fi r (t))<(l,fl r (t ))+E r (t)-E r (t ) 

(3.11) < (1, A r (to)) + ^ r (*o + - E r (t ) 

c c 

by (3.4) (with m = 0), as well as (i) and (iii) above. We now proceed by 
induction and assume that (3.10) holds on A r for some < k < \ti/l\. To 
show that the statement holds with (A; 4- 1) in place of k, we argue anal- 
ogously to (3.11). We must first show that (l,/2 r (i + (k + 1)/)) < c/2. To 
this end, we use an argument inspired in part by an idea of Grishechkin ([5], 
page 542). The idea is to consider two cases separately: the case where the 
queue length becomes zero during the interval Ik, and the case where it does 
not. If (l,p, r (s)) is never zero for s £ Ik, we can write 

rt + (k+l)l 

Sl 0+ ki,t +(k+i)i = / Vdh A r (s))) ds 

t +(k+i)i 

(l,/r(s)} -1 d*>J/c, 

t +kl 

since we have assumed that (3.10) holds for this k. By Markov's inequality 
and (ii) above, this implies that 

(% r ^ , +(1+11 „»).A r (io + kl)) < <l(i/ C)0 o),M r (<o + kl)) 
(3 12) t +ki,t +(k+i)i' i i 

<y<X,A r (*o + «)><|- 

If, on the other hand, (l,/2 r (s)) = for some s £ Ik, then all mass present 
in the system at time to + kl is gone by time s. More precisely, we have in 
this case, using m = 0, t = to + kl and /i = s — (to + kl) in (3.3), that 



(l ( 5 ro+fc; , s ,oo),A r (to + feO) = 0. 

' to+kl, s — u t +kl,t Q +(k+l)V 



t + kl,s' 

Since Sr, kLs < Sl, ku , (k , iv , we again obtain the estimate 



< 3 - 13 ) (H + ^ 0+ ( fc+1 ) ; -)' /2r(t0 + W) / =0< 4- 
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Thus, we have on A r , 



(3.14) 



< 1 '^ + ( fc + 1 )0>^(^.«o + (^)'^ + Ai0 

+ E r (t + (k + l)l)-E r (t + kl) 
c c 
^4 + 4< 

where in the first inequality, we have used (3.3) and ignored any processing 
of jobs in the second term there. The second inequality follows by (3.12), 
(3.13) and (i) above. 

Now we can complete the proof in a similar manner to that for (3.11). By 
(3.4), (3.14) and (i) above, we have on A r for any t G Ik+i, 

(1, ff(t)) < (l,fl r (to + (k + 1)1)) + E r (t) - E r (t + (k+ 1)1) 

< (l,A r (io + (k + 1)0) + E r (t + (k + 2)1) - E r (t + (k + 1)1) 



c c 
% + 4 <C - 



□ 



3.4. Weak law estimate. The following lemma provides a generalized 
functional weak law of large numbers estimate, which will be used in several 
ways throughout the paper. 

Lemma 3.4. Assume (Q.l). Let T,L>1, Z > and < e' , ?/ < 1. Sup- 
pose that g : — ► R+ is a Borel measurable function such that (g, u) < oo, 



as r 

■2- 



oo, 



(3.15) { 9 y)^(g,u) 

and such that for some p' > 0, limsup r _ >00 ((? i+p , v r ) < oo. Then 

rE r < m {-+l) 



(3.16) limsupP r sup 

r->oo \rn<\rT\ 

Proof. Define 

(3.17) c £ 



al(g,u) 



>e'\< J. 



(9,v), 
1, 



(g,f)>o, 

(g,u) = 0. 



Let Eg = e'{2acg)~ x A (1/2) and let l\ = l — e g and I2 = I + e g . Note that 
for any t G [0, L], we always have [t,t + 1] D [kie g ,kie g + Zi] for some k\ € 
{0, . . . , [L/e 9 ] }, and [t, t + Z] C [Zc 2 £g, ^2£g + Z2] for some /c 2 G {0, . . . , [L/e 9 ] }• 
For any r G 1Z, m< \ rT\ , and t G [0, L], consider the event 



(3.18) 



rE r ' m {t+l) 

E 9(v r i)-al(g,v) 

i=rE r > m (t)+l 



>£' 
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This event must be contained in the event 

rE r - m (t+l) 



1 



E g{vl)>al(g,v) +e' 



(3.19) 



i=rE r < m {t)+l 



U< 



rE r - m (t+l) 



E gK) <al{g,v) -e'\. 

k ' i=rE r - m (t)+l > 

Thus, for some k\,ki £ {0, . . . , \L/e g ~\}, the event in (3.18) is contained in 

rE^ m (k 2 e a +l 2 ) 



(3.20) 



- E g{v r i)> al 2 {g,y) + £ 

i=rE r ' m (k 2 e g )+l 



Llj- J! 5«)<a/i(ff,^) 

I j=r_B r ' m (fcie 9 )+l 



2 ' 



where we have used the definitions of s g , l\, l<i and c g . Therefore, to prove 
(3.16), it suffices to show that for j 6 {1, 2}, 



lim sup P r 



sup 

^<\rT\, k<\L/eg\ 



(3.21) 



E M. 

i=rE r > m {k£ g )+l 



alj(g,u) 



Furthermore, by summing over m < \rT\ and k < \L/e g ~\ , it suffices to show 
for j £ {1, 2} that for all sufficiently large r, all m < [rT\ and all A; < \L/e g ~\ , 



pr 



rE^ikeg+lj) 



E g(v r i) - aij{g,v) 



(3.22) 



i=rE r . m (fc£ g )+l 



< 



2(rL/ e9 l+l)(LrTj+l) 

To this end, we now fix j € {1, 2} and A: G {0, . . . , \lije g ~\ } for the remain- 
der of the proof. For better readability, we will suppress notation indicating 
dependence on k and j for certain objects, when it is clear from the defini- 
tion. 

For each r £lZ and m < \ rT\ , let 

rE r - m (ke g +lj) 



(3.23) 



1 



E 9(v r i)-alj{g, 



i=rE r ' m (ke g )+l 



^ e 
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be the event appearing in (3.22). For notational convenience, we define 
E r,m = r E r > m (ke g ) and N r > m = rE r ' m (ke g + lj) -rE r > m (ke g ). Define c+ = 1 + 
e'(8aljC g )~ l and c_ = 1 — e'(16(aljC g V 1)) , and let 

A r £ m = {lc-ralj\ < N r ' rn < [c+ralj]}. 
Let a^ 1 denotes the complement of A r ^ . For each r£?l and m < \rT\ , 

(3.24) P r (A r > m ) < P r (i^ m ) + P r (A r ' m n A r £ m ). 

We will show that each of the two terms in (3.24) is bounded above by 
rf[A(\L/e g ~\ + l)([rT] + l)]" 1 for all sufficiently large r. 

For the first term in (3.24), note that on Ajj, , either E r (r{m + ke g + lj)) — 
E r (r(m + ke g )) < [c^ralj] or E r (r(m + ke g + lj)) — E r (r(m + ke g )) > [c + ralj\ 
by definition. Inversion of the renewal process E T (-) implies that on A^ 1 , 
either 

(3.25) ^ u\ > r(m + ke g + — r(m + ke g ) = rlj 

i = E r,m + l 

or 

E r ' m +[c + ralj} 

(3.26) V] <r(m + ke g + lj) — r(m + ke g ) =rlj. 
i=E r ' m +2 

Since a r — ► a, as r — > oo [see (2.22)], and since by definition, < c_ < 1 < 
c + , (3.25) and (3.26) imply that on A r ^ m , for some C2 > and all sufficiently 
large r, either 

,3.27) £ U-± )>r, i -^i±>r* 

or 

(3-28) £ k-±)<rh- l ^^<-r« 

i=E r > m +2 ^ 

Define the random variables 



a' J a' 



(3.29) x r i =u r i t {u r < r i/ 8} and y[ = <l {u r >r .i/ 8} , 
and let 

(3.30) x:' m = £ «- Er [<])> 



E r ' m +[c^ral J \ 



(3.3i) y:> m = 2 (^ r - E % r ])> 
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E r ' m +\c + ralj\ 

(3.32) XT= £ «-E r [<]), 

i = E r, m + 2 

E r ' m +[c + ralj\ 

(3-33) y;- m = £ (yr-ETO 

i=E r > m +2 

Then using (3.27)-(3.33), we have 

P r (A E m ) < P r (|X:' m + YL' m \ > rc 2 ) + P r (|X£ m + Yl' m \ > rc 2 ) 

< V r {\X r > m \ > rc 2 /2) + P r (\X r ^ m \ > rc 2 /2) 
+ P r (\Yl' m \ > rc 2 /2) + P r (\Yl' m \ > rc 2 /2) 

< (rc 2 /2r\B r [(X r J m ) 4 ]+B r [(X r ; m ) 4 ]) 
+ (rc 2 /2y 2 (B r [(Y:> m f} +E r [(y; ,m ) 2 ]) 

< (rc 2 /2)-\[c-ral j \ 2 + ([c+raljj - l) 2 )(r 1/8 ) 4 
(3.34) + (r C2 /2)" 2 (Lc_r^j + (Lc+rnZ,-] - l))B r M) 2 ] 

+ r- 1 (c 2 /2y 2 2c + al j B r [(u^) 2 ;u r 2 >r 1 / 8 ] 

< 1 f [rr] + l / _ 1/2 2(c + ^) 2 



LrTj + lV r V (c 2 /2) 4 

The fourth inequality above follows from independence of the members of the 
sequence (resp. {y[}). The fifth inequality follows since c_ < c+, and 
the last uses Markov's inequality in the final term. Letting r — > oo in (3.34), 
we see that the quantity in the outer parentheses on the last line tends to zero 
uniformly in m [see (2.24)], and is thus smaller than r]'[A(\L/e g ~\ + 1)] _1 for 
all sufficiently large r. So ~P r (A r E m ) is bounded as desired for all sufficiently 
large r and all m < [rT\ . 

We handle the second term in (3.24) in a similar fashion. Write A r,m = 
Az U A! , where 



•Rr.tn i jyr.m , 

re 



(3.35) A r l m = l J2 9K)-ral J {g,u)<-f 

I i=E r < m +l 

and 

rE r ' m +N r ' m ,\ 

(3.36) A^ = \ giO-ral^v)^] 



2 

i=E r,m + l 
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Define the random variables 

(3.37) & r = 5«)l{ 9 (^)<rV8 } and h\ = g{v\)t {g{v r )>r i/ S] 
and let 

E r ' m +\c-ralj\ 

(3.38) G r J m = £ !■?■[#]), 

i=E r ' m +l 

(3.39) ^' m = £ W-BTO, 

i=E r ' m +l 
E r ' m +[c + ral j \ 

(3-40) G£ m = E Cflf-ETO 

i=E r > m +l 
E r ' m +\c + ral j \ 

(3.41) ^ m = E w-E r ra. 

Let C3 = e'(8c+a/j) _1 and note that by assumption, 

(3.42) K<7,f r >-<<?^>|<c 3 , 

for all sufficiently large r. So for all sufficiently large r, and all m < [rT\, 

P r (A r > m nA r E m ) 

= P r {A r ± m n A r E m ) + P r (^ m n A r E m ) 

< P r {G r l m + H r _l m + |c_mZjJ (<?, z/) - ralj(g, v) < -re'/2) 
+ P r (G r + m + il£ m + [c+raljl (g, u r ) - ralj{g, u) > re' 1 2) 

< P r (\GTl m + H r _l m \ + I [c.ral.l ({g, v r ) - (g, u)) 

+ (g,u)([c-ralj\ -ralj)\ >re'/2) 
+ P^(\GT + HT\ + \[c + ral j \({9,v r ) ~ M) 

+ (g,v)(lc + ral j \ - ralj)\ > re'/2) 

< P r (\G r _l m + H r _! m \ + c + ral jC3 + {g,v)(\c- - l\ralj + 1) > re'/2) 
(3.43) + P r '(|G^ m + H r / n \ + c + ra/ jC3 + (<?, u)\c+ - l\ralj > re'/2) 

< P r (|G r _! m + H r _! m \ > re' /A) + P r (\G r / n + H r f l \ > re' /A) 

< P r (|G r _! m | > re'/8) + P r (|G^ m | > re'/8) 

+ P r (\H r _! m \ > re'/8) + P r (|#£ m | > re'/8) 

< (re'/8)- A (E r [(G r l m f} + E r [(G^ m ) 4 ]) 
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+ (re'/8y 2 (E r [(H r J n ) 2 } + E r [(H r / n ) 2 ]) 
< (re'/Sr'dc.ral,] 2 + \_c + ral 3 \ 2 ){r 1 '^ 
+ (re'/8)- 2 ([c.ral 3 \ + [c+ral^WW,) 2 ] 

1 / [rTJ + l / _ 1/2 2(c + a^) 2 + 2c +a l j W[[h\f] 



\rT\+l\ r V' (e'/8) 4 (e'/8) 2 

In the first inequality above, we have used (3.38)-(3.41), as well as the 
fact that [c-ralj\ < N r,m < \_c + ralj\ on A r £ a . The third inequality uses 
(3.42) and the fact that c_ < 1 < c + . The fourth inequality follows from 
the definitions of C3, c_, c+ and c g , where for the first term it is necessary 
that r > VoCgje' . The seventh inequality follows from independence of the 
elements of the sequence {g\} (resp. {hi}). Using the definition of h\ and 
Markov's inequality, 

E r [(h r 1 ) 2 }=W[g(v[) 2 ;g(v[)>r 1 / 8 ] 

<r-'/V +f V>, 

which tends to zero as r — > 00 by assumption. Thus, letting r — ► 00 in (3.43), 
we see that the quantity in the outer parentheses on the last line tends to 
zero uniformly in m and is thus smaller than rf[&(\L/e g ~\ + l)]" 1 for all 
sufficiently large r. □ 

3.5. Moment estimates. We will need the following simple lemma. 

Lemma 3.5. Suppose € Mp such that as k — > 00. Suppose 

further that for some q' > and 1 < M < 00, 

(3.45) hmsup((U fc ) V( X 1+q ',Zk))<M. 

k— too 

Then for any q £ [0, q'), there exists an M q 6 [M, 00) suc/i i/tarf (x 1+<? )£) — 
M g and 

(3.46) <X 1+, ,&> — <X 1+ff ,0, k^oo. 
Furthermore, M q depends only on M , q' and q. 



PROOF. Since — > as k — > 00, (3.45) implies that for any x € 
(3-47) (t( x ,oo),0 < (1,0 = Mm (1,6) <M. 

For any rc > 1, (3.45) and Markov's inequality yield 

limsup(l (a . i0o) ,£fc) < hmsupx"( 1+9 ')(x 1+9 ',Cfc) 
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So by the Portmanteau theorem ([1], Theorem 2.1), for any x > 1, 
(3.49) {t {x ,oo),0 <limsup(l (!Ei0o) ,&> <x~^M. 

fc— >oo 

Since x 1+<3 (0) = 0> we can use integration by parts together with (3.47) and 
(3.49) to obtain 



(X 1+9 ,e>= / (l + q)xl(l {Xt00) ,£)dx 

(3.50) 



< (1 + g)M (l + J™ x-^ 1+q '-^ dx^j < 
M q = (1 + g) Aff 1 + rc- (1+<? '- g) da: 



oo. 



We can take 



(3.51) 

< oo, 

which depends only on M, q' and q. To show (3.46), first choose q\ G (g, g') 
and M qi € [M, oo) such that 

(3-52) (x 1+q \0<M qi . 

For any y > 1 that is not an atom of £, 

Kx 1+9 ,e fc ) - <x 1+9 ,OI 

+ i(x 1+9 i(,,oo),efc)-(x 1+9 i ( ,,oo),oi 

<l(x 1+9 i[o, 2/ ] > a)-(x 1+9 i[o, 2/]> OI 

+ |y-(«-«)( X 1+ « 1 l [ i j0O ),e fc )| 

+ i I/ -^-9)(x 1+ « i i [ i j0O) ,e>|. 

Since y is not an atom of £, we see that the term on the fourth line of (3.53) 
tends to zero as k — > oo. So since gi < g' and M < M 91 , (3.45), (3.52) and 

(3.53) imply 

(3.54) limsup|<x 1+ Ufc} - (x 1+q ,0\ < 2y^~^M qi . 

fc-^oo 

Since £ can have at most countably many atoms, y can be chosen arbitrarily 
large. 

□ 



The next lemma provides an upper bound for three of the moments of the 
fluid scaled state descriptors {fl r ' m (-)}- This estimate will be fundamental 
to many subsequent proofs. 
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Lemma 3.6. Assume (Q.l) and let T > 1 and < i] < 1 6e given. Then 
there exists p > and My > 1 such that for any L > 1 , 

(3.55) liminfP r (0 , :)>l-?, 

r— >oo g 

«! = ( ||(l,/2^(-))V( X ,A r ' m (-))V(x 1+P ,A r ' m (-))llL<^T}. 

[m<[rTJ J 

Proof. Choose p > such that 6p + 2p 2 < 9. By (2.31)-(2.35) and 
Corollary 3.2, using rf = n/36 there, there exists Mo > 1 such that for any 
L> 1, 

(3.56) hminfP'-(^)>l-^, 
where 

^ = (<1,//(0))V sup ||( X , / u r ' m (-))llLV(x 1+p ,/i r (0))<Mo). 
I m<|rTj J 

Define 

(3.57) Mi = 32(a Vl)M , 

(3.58) My = Mi V (M + M^l + p^ 1 )2a(x 2+2p , i/) + 4a( X 1+p , v)). 
Let L > 1 be fixed and let A\ be the event on which 

(3.59) sup \\E r ' m (- + AM )-E r ' m (-)\\ L <8aM , 

m< [rTJ 



r£ r ' m (l) 



(3.60) sup - £ X 1+P «)<2a(x 1+P ,^), 

j =rS r,m(0) + l 



m<[rT+LJ r 



rE r > m (l) 



(3.61) sup - 2 X 2+2p «)<2«(x 2+2 V>. 

i=r\E r > m (0)+l 



m<|rT+£J r 



Define A r = A\ ^ A r 2 . We will apply Lemma 3.4 three times in order to 
obtain a lower bound on the probability of A\, and hence A r '. By choice 
of p, (2.27) and Lemma 3.5 [with q' £ (p, 0) or g' G (1 + 2p, 3 + 9) there] 
imply that the assumptions of Lemma 3.4 are satisfied for each of the three 
choices g = 1, g = x 1+p and g = x 2+2p - Note that Lemma 3.4 still holds if, in 
(3.16), we replace the supremum over m < \rT\ by the supremum over m < 
[rT + L\ (see proof of Lemma 3.4). Thus, by applying Lemma 3.4 three times 
[resp., for (l,e',rj\g) equal to (4M , 4oM A (1/2), n/36, 1), (1, a(x 1+p , v) A 
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(1/2), r]/36,x 1+p ) and (1, a(x 2+2p , v) A (1/2), 77/36, x 2+2p )] and combining 
with (3.56) one obtains 

(3.62) liminfP r (^ r )>l-|. 

Thus, by (3.56) and (3.62) it suffices to show that on A r , for ali m < \rT\, 

(3.63) IKM r ' m (-)) V (X 1+p ,/2 r ' m (-))lk < Mr. 

Note that by (3.57) and (3.59), \\E r {- + 4M ) - ^(-)|| [rTJ+z < Mi/4 on 
So for c = M u l = 4M , i = and ii = [rT\ + L, conditions (i)-(iii) of 
Lemma 3.3 are satisfied on A r , which imphes that on A r , 

sup |Kl,M r ' m (-)>||L = IKl,/i r (-))llLrrj + z 

m< I rT 

(3.64) 

< Mi < M T . 

It remains to consider the (1 +p)th moment. First, consider any job i > 1 
such that C/fr" 1 < i x < i 2 for some times ii,i 2 G [0, [rTj + L]. If p, r (t*) = 
for some t* G [ti,t2], then 

X l+P (vl-Shr r - H2 ) 

(3.65) < X 1+p «-^ r -x, t *) 

= 0< x 1 +p^J 



Mi 

where the equality above follows from (3.2), by setting g = x 1+p > m = 0, 
t = \J\r~ x and h = t* - U\r- X there. If fl r {s) + for all s G [h,t 2 ], then 
on A r , 

S r tut2 >(t 2 -h) inf <p«l,/i r (s)» 

(3.66) >(*a-*i)C sup (l,/i r (s)) N ) 

Vse[ti,t 2 ] / 

^(ta-tijMf 1 , 

where the last inequality is by (3.64). So in this case, on A r , we also have 

x 1+p (<-%- V2 ) 

(3-67) <X l+p (vl-Sl lM ) 

- x I « Mi 
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Combining (3.65) and (3.67) yields that on A r , 



<x l+p 



t2-h 

Mi 



(3-68) <Ht 2 - tl) M-\oo)«)x l+p {vl) 

<(l [ (t 2 -* 1 )M 1 - 1 > oo)X- (1+P) )«)X 2+2P «) 



To complete the proof, we must show that (x 1+p >A T ' m ( s )) — Mt on A r , 
for any m' < \ rT\ and s G [0,L]. Note that since L > 1, s > 1 implies that 
p, r,0 (s) = p, r,1 (so) for some so € [0,L— 1], and fi r,1 (s) = /x r ' 2 (si) for some s\ G 
[0,L 1]. 
Thus, it suffices to consider the case m' € {0, 1} and s < 1 and the case 
m' G {2, . . . , [rTj } and s G [0, L] . If m' G {0, 1} and s < 1, we have on A r , 

(x 1+p ,A r ' m 'W) 

= (x 1+p ,A r K + s )) 

= (x i+p ( ._^ (m ' + s) ),^ (0)) 

j r£ r (-m'+s) 

(3.69) r-E r (m'+s) 

<(x 1+p ,A r (0)) + - £ x 1+p (v r i) 

r i=l 

< Mo + £ J £ x 1+p «) 

m=0 V i=rE r > m (0)+l ' 

<M + 4a(x 1+? » 

< Mr, 

where the second equality above uses (3.2), with g = x 1+p , m = 0, t = 
and h = m! + s there, the first inequality follows by simply ignoring any 
processing which has taken place by time m' + s, the second inequality 
follows from the definition of A r for the first term and by including additional 
summands in the second term, and the third inequality is by (3.60) and the 
fact that m' < 1. 
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Last, we consider the case m' G {2, . . . , [tT^J } and s G [0, L\. Let N = 
[m! + s\ - 2 and note that N G {0, . . . , \ rT + LJ - 2}. Note also that by 
replacing m by [m' + s\ — m — 1, the following estimate holds: 

" ( Mi n 1+p 



o V L m ' + s J ~~ 771 — i 



m=l 

(3.70) 



m J 



<Ml +p (l + J™ x-( 1+ ri dx 

= Ml +p (l+p- 1 ). 
For m' G {2, . . . j [rT^J } and s G [0, L] we have on A r , 
(X 1+p ,//' m '(s)> 



rE r (m'+s) 

r i=i 

N /, rE r < m {\) \ 

< M o + E - E x 1+p « - 

m=0\ j=rJ5 r . m (0)+l / 
r^ r (m'+s) 

(3-71) +- E X 1+p «) 

i=rS r (AT+l)+l 

* // Ml \(1+ P )i \ 

^ \\[m' + s\ -m-lj r . = ^ (Q)+1 ^ 

+ E ; E 

< M + M 1 1+P (l + p- 1 )2a( X 2+2p » + 4a( X 1+p , 

< Mr, 

where the first inequality is obtained as in (3.69), the second follows by 
ignoring any processing of jobs in the third term; the third inequality is by 
(3.68) (with t\ = m+ 1 and t<i = m' + s) for the second term and by including 
additional summands in the third term; the fourth inequality is by (3.60), 
(3.61) and (3.70) and the last is by (3.58). □ 
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3.6. Combining estimates. To aid the reader, we now provide a lemma 
that combines several of the preliminary results obtained in Sections 3.2-3.5, 
such that they hold simultaneously, with respect to specific choices of various 
constants. It is analogous to Lemma 5.2 in [6], in that subsequent references 
to results obtained thus far will predominantly reference this lemma. We 
will need one more definition. 

Definition 3.7. Let C = {g G Cj(R+) : g(0) = g'(0) = 0}, and let C = 
{g GC:g has compact support}. Let V = {g% }keN C C be a countable subset 
such that for any g EC, there is a sequence {gk}kLi C V that together with 
{g^^Li is uniformly bounded and satisfies 

(3.72) g k — >g and g' k — > g' , k->co, 

pointwise on R + . 

An example of how such a set V CC can be chosen is given in the proof of 
property (iii) for Theorem 5.1 in [6]. Recall that M„ is defined in Definition 
2.1 and the metric d[-, •] is defined in (1.4). 

Lemma 3.8. Assume (Q.l) and let T, L > 1 and < rj < 1 be given. Let 
{e n }%Li C (0,1) be a sequence such that e n J, 0, as n — ► oo. Then for each 
n G N, t/iere exisi strictly positive constants l n , p, Mt, K n , j n , N n , r n and 
events {B^r^n, such that: 

(i) P r (B r n )>l-7], for allr>r n , 

(ii) r n+ i > r n and B r n+1 C fi^, and 

(iii) /or eac/i r G 7£ and eac/i m < [rTj , on the following hold: 



(3.73) 




(3.74) 



sup E r ' m {t + Z n ) - E r ' m (t) < — , 



fG[0,L] 



(3.75) 



sup 

k<n 

i=i,...,Li/i n J 



r 



1 



rE r '™(jl n ) 

9k{v r i)-ajln{gkiv) <e, 



i=rE r ' m (Q)+l 



for g% G V, 

(3.76) sup (l,^ m (t))V(x,/2 r ' m (t))V( X 1+p ,/i r ' m (t))<M T , 



te[o,L] 



(3.77) inf d[//(0U]V|(x,//(0)>-(x,£)l<£. 



(3.78) 



71 



2Mt ' 



DIFFUSION APPROXIMATION FOR A PROCESSOR SHARING QUEUE 31 



(3-79) sup (1 ^(o)) < ^, 

(3.80) 7n<^A^, 

(3-81) sup |(x,/i r ' m (<)>-(x > /i r ' m (0))|<^ ) 

(3.82) N n = \64M%/(K n e 2 n )] + 1, 

rE r > m (t+l„) 



(3.83) 



SU P - X! %K„,(fc+l>n)(<) 

te[o,L] r 



i=r£ r ' m (t)+l 

t 7, 



- ~jf ( 1 [(fc-l/2)K„,(fc+3/2)K n )^) 

/or all k G {0, 1, . . . ,JV n }. 

Proof. Let fj = i]/6. For each n G N, choose < l n < e n j (32a) such that 
ln+i < In and let be the event on which both (3.74) and (3.75) hold. 

Let Bn,l = n m <|rTj ^n'T and let = Hn<n The first § oal is to show 

that for each n G N, 

liminfP r (^ 1 ) > 1 -77. 

r — >oo 

Fix n G N and consider any h<n. Since (3.74) (with n = h) holds on B 7 h l 

for each m < \rT\, and since (3.75) (with n = h) holds on B r h 1 for each 
m < \rT\ , we must apply Lemma 3.4 1 + h [L/l n \ times in order to estimate 
the asymptotic probability of B r h l . More precisely, we must apply it once to 
guarantee (3.74) and also once for each case k < h and j < \L/l n \ in (3.75). 
Observe that by (2.25), the assumptions of Lemma 3.4 are satisfied with 
/ = In, e' = £n/16, rf = 77(71(1 + n[L/l n \))~ 1 , p' > arbitrary and g = l. Also 
observe that for each k = 1, . . . , h and j = 1, . . . , [L/l n \ , the assumptions 
of Lemma 3.4 are satisfied with l=jl n , = £h, v' = v( n (^ + n \_L / l n \))~ l , 
p' > arbitrary and g = g^.. Thus, 1 + h[L/l n \ applications of Lemma 3.4 
yield 

(3.84) Itataf P-^) >!-(! + WUJ) ^,,^, > 1 - * ■ 
This implies that for each n £ N, 

(3.85) liminf V r (B r n l ) > 1 - J2 ~ = 1 ~ »7- 

n=l ra 

Note that by definition, 1 C JB£ 1 for n G N. 
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Next, let p and Mt be the constants given by Lemma 3.6 and let B^'™ 1 
be the event appearing in (3.76). Let B\ = r\ m <\rT\ B^ m . Then by Lemma 
3.6, 

(3.86) limmf V r {B r 2 ) >l-fj. 

For each n G N, let B T n 3 be the event appearing in (3.77). Recall that 
"~" denotes equivalence in distribution. By (2.35) and the Skorohod repre- 
sentation theorem, there exist random pairs (@ r ,X r ), (G,X), taking values 
in A4p x and defined on some common probability space (f2, J^, P), such 
that (e,X)~(0,(x,e» and (9 r ,X r ) ~ (/2 r (0),(x,/I r (0))) for each reU 
and such that as r — ► oo, 

(3.87) (6 r ,X r ) -^(6,X) a.s. 
So by (2.30) and (3.87), 

liminfPVB;,) 

(3.88) =liminfP( inf d[G r ,£] V \X r - ( X ,0\ < ^ 

> 1 ~fj. 

Since {e n } is a decreasing sequence, B^ +1 3 C B^ 3 for n G N. 

Next, for each n G N, we use (2.34) (with e = e n /4) and (2.31) to choose 
K n > and -M^, > so that (3.78) holds, K n +i < K n for n G N and 

(3.89) P(jup (l [X)X+2Kn] ,e) V (l [Mn ,oo),0) < ^) > 1 " £f 

For k G {0, ... , [M n /K n ] - 1}, define I k = [kn n , {k + 2)n n \ and let I\M n / Kn \ = 
[M n ,oo). Note that for any x G R+, + ft n ] C Ik for some fc €E {0, 
[M n /K n ]}. Consider the random measures G r ,0 defined above such that 
(3.87) holds. Since for each n G N, {Ik} is a finite set of closed intervals, 
a trivial generalization of the Portmanteau theorem ([1], Theorem 2.1) to 
finite measures yields for each n G N, 

(3.90) limsup max (l/ fe ,O r )< max (l/ fc ,0) a.s. 

rwoo A;<[Af n /re n ] fc<[Af„/K„] 

Since 6 ~ 6, (3.89) and (3.90) imply that for each n G N, 



(3.91) >liminfP( max (l A ,e r )< — 
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Let B r n 4 be the event appearing in (3.79). We have by (3.91), 
liminfP r (5; 4 ) 



(3.92) = liminfP sup {t [x , x+Kn] , @ r ) < ^ 



7] 

~ ~ 2^2' 



Letting B r n 4 = C\ i<n -B[ 4 , we see that for each n G N 

2i 2 



(3.93) liminf P^, 4 ) > 1 - £ J* > 1 - 



i=l 

By definition, -B^ +1 4 C B T n 4 for n G N. 

For each n G N, choose 7 n > such that (3.80) holds and such that 7 n +i < 
7„ for n G N. Let -B^'™ be the event appearing in (3.81) and let B r n 5 = 
rim<|rT| ^n'T- Then by Corollary 3.2, using = fj and 7 = 7 n , we have for 
each n G N, 

(3.94) liminf P r (B£ )5 ) > 1 - fj. 

Since {7 rt } is a decreasing sequence, B n+1 5 C .B^ 5 for n G N. 

Finally, define N n as in (3.82) for each n G N, and for A; G {0, . . . , iV n }, 
choose G C&(]R + ) such that g% is nonnegative and such that for all x G R+, 

l[fc« n ,(fe+i)«„)(^) < 5^0*0 < l[(fc-i/2) Kn ,(fc+3/2) Kn )(2;)- Note that for fc = and 

X GM+, l[(fe-l/2)«„,(fe+3/2)K„)(^) = l[O,3/2«„)O0- Let ^nj^ be the event in 

(3.83) and let B r nfi = f] k <N n ^m<\rT\ B%fi k - Define 

(3.95) G; 6 =n sup sup- J2 9U<)<fi9%,v) ■ 

fc =o lm<LrTj te[o,L] r . =r ^ rjm(f)+1 » J 

Fix n G N and consider n<n. Observe that by (2.25), we have that for each 
k G {0, . . . ,N n }, the assumptions of Lemma 3.4 are satisfied with I = l n , 
e' = {e n /lQ){g^,v), = rj(n(N n + l)) -1 , p' > arbitrary, and g = g^. Thus 
N n + 1 applications of Lemma 3.4 yield 

(3.96) liminf P r (G r ne ) > 1 - ^ + ]\ > 1- 1 

y ' r-*oo K n ^> ~ n(N n + l)~ n 

where we have used the fact that l n < e n /(32a), and the second inequality 
follows from (3.82) since {n n } and {e n } are decreasing sequences. Since 
G r n 6 cBJ 6 for each h < n, we have 

(3.97) limmfP^ i6 )>l-^. 
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Letting B r n>6 = f| 

?i<n^n,6' ®^ that B n ^ Q C B n Q for ti £ N, and 

n 

(3.98) hminf P^ 6 ) > 1 - £ J = 1 - rj. 

n=l 

For each r £ ft, define 5; = BJ^ n J3J D B r n Z n B^ >4 n ££ 5 n J3^ 6 and note 
that by construction, B r n+1 C for n € N. Then by choice of 77, we see that 
for each n £ N, there exists r n £ ft such that r n+ i > r„ and r > r n implies 

(3.99) P r (B r n )>l- V . □ 

3.7. Fluid scale analysis. In this section we build on the results of Sec- 
tions 3.1-3.6 to prove several important properties of the fluid scaled state 
descriptors {/J r,m (-)}- These properties will be used at the end of the sec- 
tion to prove the desired relative compactness property (see Corollary 3.16). 
Lemmas 3.9 and 3.10 give an upper bound for the amount of mass that 
p, r ' m (t) can have concentrated near zero. They are analogous to Lemma 5.5 
in [6]. In Lemmas 3.11 and 3.12 we give an upper bound for the shifted fluid 
scaled queue length {1, p, r,m (-)) when the shifted fluid scaled initial workload 
(x,M r,m (0)) is small, and a lower bound for (l,p, r > m (-)) when (Xi/J r ' m (0)) is 
bounded away from zero. These two lemmas are analogous to Lemmas 5.3 
and 5.4 in [6]. We then combine the four aforementioned lemmas to prove 
Theorem 3.14 which provides an oscillation bound for the sample paths of 
the processes {/J r ' m (-)}> and yields Corollary 3.16 as a consequence. We begin 
with the following technical lemma. 

Lemma 3.9. Assume (Q.l) and let T, L > 1 and < rj < 1 be given. Let 
{e n }^ =1 C (0, 1) be a sequence such that e n { 0, as oo. Let l n , p, Mt, 
k„, 7„, N n , r n be the constants, and {B^} be the events, given by Lemma 
3.8. Fix n £ N, r £ ft, m < [rT\, and t £ [m,m + L]. Define the following 
random times: 

r £n =sup| S £[0,t]:(l, / u r ( S ))<| 

f 8M£ 
t = max 0,r £n ,t 



where we define T £n = — oo for the supremum of the empty set. Then on B r n , 
(3-100) (1(0iKn]( ._^ ) ,^ (r)) <^. 

Proof. We only consider realizations in B r n in the following. We treat 
each possibility for r separately, and suppose first that r = 0. Then by (3.79) 
we have 

(3.101) (l (0jKn] (- - K,t)^ r (r)) < sup (l [x , x+Kn] ,(i r (0)) < ^. 
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Next suppose r = r £n > 0. Then there exists a r' € [(r En — l n ) V 0, t £ J such 
that (l,/2 r (r')) <£ n /8. This implies by (3.4) (with m = 0) and (3.74) that 

(l (0A?i] (.-S;,),^(r))<(l,^(r)} 

(3.102) 

<^ + E r (r' + l n )-E r (r') 

o 

- £n £n 

_ i ie ? 

Lastly, suppose T = t- (8M T /e n ). Let t" = t- (4M|/e n ). Then since r" > 
r £n , we have for s € [r",i] that (l,/2 r (s)) > e n /8. This implies by (3.76) that 

r l -, i-r" 4M T 

(3.103) S r T ,>S r T „,= (l,fl^ s ))-Us>— — >— 

So by (3.103), Markov's inequality and (3.76), 

(l (0 , Kn] (• - K,tW^)) ^ <%oo)(- - S r T>t ), fi r (r)) 

< (^{4,M T /e n ,oo),V r ( T )) 

(3.104) 



- T 



□ 



The next lemma gives, on B^, an upper bound for the amount of mass 
that p, r,m (t) can have concentrated near zero for t £ [0, L]. 

Lemma 3.10. Assume (Q.l) and let T, L > 1 and < i] < 1 6e given. Let 
{e n }^ = i C (0,1) 6e a sequence such that e n [0, as n — > oo. Let l n , p, Mt, 
k„, 7„, iV n , r n be the constants, and {B T n } be the events, given by Lemma 
3.8. Then for each n £ N, rGK and each m < [rT\ , on B r n , 

(3-105) sup (l [0)l(n ],A r ' m (<)><v- 

Proof. In this proof we only consider realizations in B r n . Fix n £ N, 
r E 7£, m< [rT\ and i € [m, m + L], and let r be defined as in Lemma 3.9. 
We must show that 

(3-106) <%*j,M r (t)><y. 

If t = t, then the result follows from Lemma 3.9, since fl r ' m (t) does not 
charge {0}, and S T T T = 0. Thus, it suffices to consider the case r < t. Consider 
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an arrival i such that U[ r -1 G (t, i] and v\ > (64M|,/e^) + 2ft n . Since t — T< 
8M|./e n , and since (l,j2 r (s)) >e n /8 for s G (7", i] , we have 

(3.107) S r UTr - ljt < S r T>t = /WoO)" 1 ds<(t-r)^< 

Thus, u[ - t > (64M|/4) + 2k„ - (64M|/4) > Kn . Recalling that 

jV n = [64M|/(«; n 4)] + i 5 this implies that for > iV n + 1, 

(3.108) l [fcKni(fc+1)Kn) «)l (0iKn] « - ^j. r _i jt ) = 0. 

Next, consider two jobs i < j for which t>[,i>J G [fc«; n , (A: + 1)k„) for some 
integer fe G {0,...,N n }. If C/Tr -1 , C/jr" 1 G (r,t] and C/Jr" 1 - f/fr" 1 > Z n , 
then, using the fact that (l,ff(s)) > for s G (r, i] again, we have at time t 
that 

(uj - S UVr -i t ) - (u[ - S UVr -i t ) 

J % 

= '5{/ r 'r- 1 ,C/ r r- 1 + U i ~~ U i 

> 1 



sup se [ , i ](l,^ r (s)} 

where the last two inequalities are by (3.76) and (3.78), respectively. This 
implies that at most one of 

is nonzero. So for each fe G {0, . . . , N n }, all jobs z satisfying U[r G (r, i], 
u£ G [kn n , (k + 1)k„) and v\ — S 1 U r r _ 1 ( G (0, /c n ] must also satisfy C/J'r^ 1 G 

(s, s + Z n ], for some s G [r, (i — Z„) V r]. (Note: s is random in general.) This 
yields the following estimate at time t, for each k G {0, . . . , N n }: 

1 rE-(t) 



Y 1 [fe«„,(fc+i)«„)( w i) 1 (o,K„](' y i - S'f/r 

r£ r (s+Z n ) 



i=r£ r (T)+l 



(3.109) < sup - ]T %K„,(im)K„)(<)%o,M(< - § urr-\t) 

«6[T,(t-I„)VT] i = r ^ r(s) + 1 

< SUP SUP - l[fc«„,(fc+l)«n)(«i)- 
m<[rTJ se[0,L] r i=r £r,m( s ) +1 
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Using (3.2) (with m = 0), we have that on B r n , 
(t[ , Kn ],(f(t)) 

= (i ( o iKn] (- - K,tU r ^)) + ~ E H°,k«M - 

i=rE r (r)+l 

£ oo 1 rE-(t) 

k=0 i= r E r (T)+l 

N n 1 rE r - m (s+l„) 

<4 + E SU P SU P ~ E ^[kK n ,(k + l)K n )K) 

4 k=0 ™<lrT] se[0,L] r i=r Er, m{s)+l 

^ -J + H -7r(M(k-l/2)K n ,(k+3/2)K n ),v) 
k=0 

I oo \ 

^ -j + y^z^%fc-V2)^,(fc+3/2) Kn ),iy 
~T + T' 

where the first inequality is by Lemma 3.9, the second inequality is by 
(3.108) and (3.109), the third inequality is by (3.83) and the last line follows 
since v is a probability measure. □ 

The next lemma gives an upper bound on [0, L] for the process (1, fL r ' m (-)), 
on the event that (x>A r,m (0)} is below a threshold. 

Lemma 3.11. Assume (Q.l) and let T, L > 1 and < i] < 1 be given. Let 
{e n }%Li C (0, 1) be a sequence such that e n { 0, as n — > oo. Let l n , p, Mt, K n , 
In, N n , r n be the constants, and {B^} be the events, given by Lemma 3.8. 
For each n G N, r G K and m < [rT\ , define the event DHj™ = {{x, /i r ' m (0)} < 
7 n /2}. Then for each n G N, r G TZ and m < [rT\ , on B r n n Dlj™ , 

(3.110) sup (l,/2 r ' m (t))<2e n . 

te[o,L] 

PROOF. Fix n G N, r G 7£, and to < [rTj . We must show that on B r n n 

(3.111) sup (l,p, r (t)) <2e n . 

t^[m,m+L] 

Let A r = B r n fl DJj™, and let to = m, t\ = L, c = 2e n and / = l n . Then con- 
dition (i) of Lemma 3.3 holds by (3.74) and condition (ii) holds by (3.80), 
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(3.81) and the definition of Dlj™. To see that condition (hi) holds, observe 
that on Bllr\D r ' m , 

= (l[o, K „],^ m (0)) + (l( Kni0 o),^ m (0)) 
(3-112) <^ + _L( X) ^ m(0)) 

< £n + J_2n 

where the first inequality is by Lemma 3.10 and Markov's inequality, the 
second uses the definition of D'Hf™ and the last is by (3.80). So the result 
follows by Lemma 3.3. □ 

The next lemma provides for each m < \rT\ , a lower bound for the process 
(1, /2 r ' m (-)) on the event where (x, /2 r,m (0)) is above the threshold 7 n /2. 
A consequence of this is an upper bound, on this event, for the rate at 
which S r t '^_ h can increase as a function of h. 

Lemma 3.12. Assume (Q.l) and let T, L > 1 and < i] < 1 6e given. Let 
{e n }%Li C (0, 1) 6e a sequence such that e n j 0, as n — > oo. Let Z n , p, My, K n; 
7 n , iV n , r n be the constants, and {B^} be the events, given by Lemma 3.8. 
For each n£N, r € TZ and m < \rT\ , let Dlj™ be the complement of D?j™, 
that is, i)Z>™ = {(x,P- r ' m (0)) > 7n./2}. Then for each n € N, there exists T n > 
such that for all r £ 7£ and m < |_r-^J , on B T n n -D^ 

(3.113) inf <M r ' TO (t)>>^- 

*G[0,L] 1„ 

and 

(3.114) sup Slf +h <hY n 

te[0,L-h] 

for any < h < L. 



•\r,m 



B r n D r,m 
n 1 1 7n 



Proof. Fix n G N, r GlZ and m < We have for any if > 0, on 

te[o,L] 1 J v . / 
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< inf (K(t [0iK] ,fi^(t)) +K-P( X 1 +P,fl^(t))) 
te[o,L] 1 J 

< inf K(l,p, r > m (t)) + K~ P M T} 
te[o,L] 

where the first inequality is by (3.81) and the definition of Dy™, and the 
third is by (3.76). So for any K > 0, on B r n n fy™, 

(3.115) WW0)^S/tp*. 

te[o,L] K 

Setting K = K n large enough so that K~ p Mt < 7n/5, and letting T n = 
K n {^f - If)- 1 proves (3.113). To prove (3.114), we have by (3.113) on B' n n 

f)r,m 

for < h< L, 

sup S[f +h <h sup y{{lrf>™(t))) 
te[o,L-h] te[o,L] 



<h[ inf (i,pr> m {t))) 1 <hT n . 

te[o,L] 



□ 



We are now ready to apply the previous four lemmas to prove an oscilla- 
tion bound on B T n for the sample paths of /2 r,m (-), m < \rT\ . We use Lemma 
3.11 to prove it on B r n r\Dy>™ and use Lemmas 3.10 and 3.12 to prove it on 

B^CiDlj™. Recall that T)l{Mf) is the Skorohod space of r.c.1.1. functions 
defined on [0,L] and taking values in A4p. In order to state the theorem, 
we define a modulus of continuity on ~Dl(A4f), with respect to the metric 
d[-, •] as follows. 

Definition 3.13. For any L > 1, £(•) S ~D l (My) and 5 > 0, define 

(3.116) w L (C(-),<5)= sup sup d[t(t + h)X(t)}. 

te[o,L-8] he[o,s] 

Theorem 3.14 (Oscillation bound). Assume (Q.l) and let T, L > 1 and 
< 77 < 1 6e given. Let {e n }^ =1 C (0, 1) 6e a sequence such that e n J, 0, as 
n — » 00. £e£ Z n , p, Mt, K n > In, N n , r n be the constants, and {B^\ be the 
events, given by Lemma 3.8. Then for each e > there exists 6 > and 
n £ 6N such that for all n> n £ and all r € TZ, on B r n , 

(3.117) sup w L (ft r ' m (-),5)<e. 

m< \rT\ 

PROOF. Fix e > 0. Choose fc £ eN large enough so that J2T=k e +i 2_fc < 
e/4, and let £T e = (max^ =1 ||( 5 f )'|| 00 V 2). Recall that Q = {g G k }f =1 U {/ijftgLi 
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is the set of functions used to define the metric d[-, ■] [see (1.4)]. Let Q e = 
{0fc}*LiU{/i£}£Li so that for any / e a £ , ll/IU < 1 and H/'lloo < i^ E . Choose 
n £ £ N large enough so that e„ e < e/(8k e ), and let 

(3.118) 5 = min(-,Z ne , 



U' " e '8r ne m t k £ (K £ viyr ne 

Fix n > n e , r € 1Z and m < I rT I , and consider the event DL' m . On BT. D D^' m , 
we have for all / G £/ £ , 

sup sup K/,/2 r ' m (t + ^))-</,A r * m (<))l 
te[0,L-5] /iG[0,5] 

(3.119) <2 sup |(/,A r ' m W)|<2||/||oo sup <l,A r,m (t)) 

te[o,L] te[o,L] 



The third inequality above follows by Lemma 3.11, since B r n n DJj™ C ££ e D 
-D?^ 71 for n > n £ [see Lemma 3.8(h)], and by the fact that ||/||oo < 1 by 
definition. The last inequality above follows by choice of n £ . We must show 
that the above estimate also holds on BZ. n D r > m for all f (zG £ . First, observe 

that on B r n n D^'™, a first-order Taylor expansion of / gives the following 
estimate for all < h < L, t € [0, L — h] and y S (i?j' t , h , oo): 

(3.120) \f(y - S r t f +h ) - f(y)\ = | - S r t £ h f'(w y )\ < fcr„J/'||oo < hT n K £ , 

for some w y € [y — ££'£? ft ,y], where the first inequality follows by Lemma 

3.12, since B T n n DZj™ C B r ns n -D^" 1 for n> n £ , and the second inequality fol- 
lows by definition of Q e . Now subtracting (/, fl r ' m (t)) from both sides of equa- 
tion (3.2) and using the fact that (l (0lOO )/)(- - = 1 (^ h .oo)(0/(' ~ 
<^t+/J for i G [0, L - /i] yields that on B r n n 

K/,/i r ' m (t + M>-</^ r ' ro (t))l 

-(l mr+h] f,fi r > m (t)) 

rE r < m {t+h) 

(3-121) +- £ (l(o,oo)/)K-^-- 

i=r\E r > m (i)+l 



< (|i(5[- h ,oc)0)(/(- " §& h ) ~ f(-))lfi r ' m (t)) 

+ ||/||oo<%hr^i, ^W) + \\f\UE^(t + h)- E r ' m (t)) 
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< hT ne K £ (l,^ m (t)) + <l [0)MW ],/i r ' m (f)> 
+ (E r > m (t + h)-E r > m (t)), 

where the first inequality is by Lemma 3.12 (again using the fact that B r n n 
i)^™ C BL^nDJj™ for n > n e ) and the second inequality is by (3.120) and the 
fact that ll/Hoo < 1- Taking the supremum over h £ [0,(5] and £ £ [0, L — <5], 
we see that on BL n .&' m , 

sup sup \(f,ii r > m (t + h))-(f,ji r > m (t))\ 

t£[0,L-8\ he[0,6] 



< sup [6T ne K £ (l,fi r > m {t)) 
te[o,L-6] v 

(3.122) + (l [0t5Fne] ,fi r >™(t)) + (^ m (i + 5) - 



te[0,L-<5] \ oM T k £ 



+ (t [0 , Kne] ,fi r ' m (t)) + + U - ^ m (0) 



< — - — Mt + ^ + ^ 
- 8M T k £ 2 16 

e e e 



\Qk £ 12%k £ 

e 

where the second inequality is by (3.118) and the third is by (3.76), Lemma 
3.10, and (3.74) (once again using the fact that B r n n D~L>™ C B r Us n D~L>™ 
for n > n £ ). So the desired estimate (3.119) holds on both B r n n Dlj™ and 

i?£ Pi -D?^ 71 and therefore on B r n for all n> n £ , r £lZ and m < [rT\ . Thus, 
combining (3.119) and (3.122) with the definitions (3.116), (1.4) and the 
definition of Q e , we have on B^, for all n> n £ and r€K, 

sup w L (//> m (-),(5) 

m< [rTJ 

= sup sup sup d[(i r > m (t + h),ji r ' m (t)] 
m<|rTJ te[0,L-<5] he[0,<5] 

< sup sup sup (^2-\\{gln r ' m {t + h))-{glff' m {t))\^) 
m<[rT\te[0,L-S] he[0,S]\ k=1 



]T + sup \{h%,fi r > m (t + h)) 



k=k e + l fc GN 
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< — - 



2 2L 



e 

+ - 



2k F 



(hlfi r ' m (t))\ 



□ 



The preceding theorem has the following important consequence, which 
will be used to identify certain measure valued paths which approximate the 
sample paths of the fluid scaled state descriptors {fi r ' m (-)} on the events B r n . 
This will be the first important step in our state space collapse argument. 
We will need the following definition. 

Definition 3.15. Assume (Q.l) and let T,L>1 and < r) < 1 be given. 
Let {e n }^ =1 C (0, 1) be a sequence such that e n [ 0, as n — > oo. Let l n , p, Mt, 
Kni In, N n , r n be the constants, and {B^} be the events given by Lemma 
3.8. Then for each r G H, we define 



J sup{n G N : r n < r}, r>r±, 
\l, r<r 1} 



(3.123) n(r) 

and we define 

^£ = {C(0 e D L (M F ) : C(-) = //' m (0M on [0, L], 

(3.124) 

for some uj G -B^( r ) and m < \rT\ } . 

Note that since r„ — ► oo as n — ► oo (see Lemma 3.8), n(r) is finite for each 
r G 1Z. Note also that for r > n, ^£ is nonempty since by Lemma 3.8, 
F-(2% r) )>l-»7. 

Corollary 3.16 (Relative compactness). Assume (Q.l) and Zei T,L> 
1 and < rj < 1 &e given. Let {e n }^ =Q C (0, 1) be a sequence such that e n J, 0, 
as n — > oo. Let IZClZ be a subsequence and suppose {C r (')}fe7l c D/^A^f) 
zs a sequence of paths such that for each f > ri, £ r (-) G T/ien {C r (")}fe7?. 
is relatively compact in D^(.Mf). Moreover, any limit point £(-) o/ i/ie se- 
quence is continuous. 



Proof. Define 

c T = {(eM F :(i,0v(x,0<M T }. 

Since {x, C) < Mr implies (l[A',oo)i C) < Mt/K for any if > 0, we have 

SU P (%,oo)>0 ->0 asAT^oo, 
CeC T 
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which implies that Ct C A4f is relatively compact ([10], Theorem A7.5). 
Let Ct be the closure of Ct and observe that by Definition 3.15 and (3.76), 
C r (t) £ Ct for all r > r% and t £ [0, L\. Since n(f) — > oo as f — > oo, Theorem 
3.14 implies that 

(3.125) lim lim sup w L (C ? (•),£) = 0. 

Thus, the relative compactness of {C r (")}re7l follows from [4], Chapter 3, 
Theorem 6.3, by noting that the modulus of continuity used there is bounded 
above by w^(-, •), and that the result stated there still holds if one replaces 
Doo(A^f) by D L (At F ) and T > there by L. By (3.125) and the definition of 
w l('j ■)> we see that any limit point £(•) of the sequence must be continuous. 
□ 

4. Diffusion limit. Having established the desired relative compactness 
property for sample paths of the processes {ji r ' rn (-)}, we now turn to the 
second main task in our strategy for proving Theorem 2.3. In Section 4.1, 
we identify limit points of the sequences specified in Corollary 3.16 as fluid 
model solutions on [0,L] (see Definition 4.2 and Lemma 4.3). We then show 
that for large r, sample paths of the processes {fi r ' m (-)}, can be uniformly 
approximated by fluid model solutions which are arbitrarily close to being 
in steady state (see Lemma 4.4 and Proposition 4.5). This leads us to the 
proof of state space collapse and finally to the proof of Theorem 2.3, both 
of which are contained in Section 4.2. 

4.1. Local fluid approximations. 

Definition 4.1. Assume (Q.l) and let T, L > 1 and < ?/ < 1 be given. 
Let {e n }^ =1 C (0, 1) be a sequence such that e n [ 0, as n — > oo. Define 3§l to 
be the set of all ((■) £ D^(A / Jf) such that there exists a subsequence TZclZ 
and a sequence {C r (')}fe7^ c D^A-^f) satisfying £ r (') € ^£ for each r £1Z, 
and 

(4.1) CT)^C(-) asf-oo. 

Definition 4.2. Given 1 < L < oo, a fluid model solution on [0,L] for 
the critical data (a, v) is a function £(•) : [0, L] — > Aip such that: 

(i) ((■) is continuous, 

(ii) (1 {0} ,C(*))=0 for allte [0,L], 

(iii) if C(0) + 0, then ((t) ^ for all t £ [Q,L] and £(■) satisfies 

(4.2) (9,m) = (gX(0))- f t{ -^^-d S + at(g,u), 

Jo (l,C(s)> 

for all t £ [0, L] and all g £ C = {g £ Cj(M+) : 5(0) = 0, g'(0) = 0} and 
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(iv) if C(0) = 0, then ((t) = for all t G [0,L]. 

We note that Definition 4.2 differs slightly from the definition of a fluid 
model solution for critical data (a, u) given in [6], Section 3.1. Besides the 
fact that we only consider fluid model solutions defined over finite time 
intervals here, the time t* used in [6] (the first time at which a fluid model 
solution reaches the zero measure) is not present in our definition. In fact, 
it is shown in [6] that if £(0) ^ 0, then t* = oo. Indeed, if we replace [0, L] 
by [0, oo) in Definition 4.2 above, we obtain an equivalent definition to that 
given in [6]. Similarly, by restricting a fluid model solution as defined in [6] 
to the finite interval [0,L], one obtains a fluid model solution on [0, L\ as 
defined here. The next lemma asserts that the elements of £$l are in fact 
fluid model solutions on [0, L] . 

Lemma 4.3. Assume (Q.l) and let T > 1 and < rj < 1 be given. Let 
{e n }^ = i C (0, 1) be a sequence such that e n 10, as oo. Then there exists 
q € (0,p) and 1 < M q < oo such that for any L > 1 and any £(■) £ 

(i) su Pte[0 , L] (i, c(t)) v ( x , m v ( X l+q , at)) < M q , 

(ii) (X)C(')) * s constant on [0,L] and 

(hi) a~) * s a fluid model solution on [0, L\ for the critical data (a,u). 
Proof. Fix L > 1 and let l n , p, M T , 

K n> Tru -^?i5 T n be the constants and 
{B' n } be the events given by Lemma 3.8. Fix £(•) £ 3Sl and let 1Z C 1Z be 
a subsequence and {( r (•)}- g ^ C D^(AIf) be a sequence such that £ r (-) G ^£ 

for each r ElZ and C( - ) — 1 -> C(-) as f — > oo. Recall that for each r eH, C r ( - ) 
is a realization of fL r ' m {-) on 5^-^ for some m < \ fT\ . This fact will be used 

throughout the proof. By Corollary 3.16, ((■) is continuous, so as f — > oo, 
we have 

(4-3) ||d[CT),C(-)]llL-o. 

By (3.76), 

(4.4) limsup sup (l,C(t)) V ( X ,C(t)) V ( X 1+p ,C(t)) < M T . 

r-Kx> te[0,L] 

Combined with (4.3), this implies that 

(4.5) sup <l,C(i)>< Mr- 

te[o,L] 

Similarly, combining (4.3) and (4.4) with two applications of Lemma 3.5 
(with q' = p and M = Mt) we see that there exists q G (0,p) and M q € 
[Mr, oo) such that 

(4.6) sup ( X ,at)) V( X 1+q ,at))<M q , 
te[o,L] 
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and such that for each t G [0, L] , 

(4.7) |<X,r(t)) -<X,C(t)>|-0 asf^oo. 

Together, (4.5) and (4.6) imply (i) above. Notice that the constant M q given 
for p and Mt by Lemma 3.5 does not depend on L, and that M q > My > 1 
[see (3.51) and Lemma 3.6]. 

Next, we have for any tG [0, L], 

l(x,C(*))-(x,C(0))| 

<liminf(|( X ,C(t))-(x,C r "W)l 

(4-8) + \(x,C(t)) - ( X ,C(0))\ + \(x,Cm - (x,C(o))|) 

< lim inf n ^ 

f^oo 4 

= 0, 

where the second inequality is by (4.7) for the first and third terms and 
(3.81) for the second term. The last line follows since n(f) — > oo, as f — > oo, 
and 7 n — > 0, as n — > oo (see Lemma 3.8). Thus (x,C( - )) i s constant on [0, L], 
which proves (ii) above. 

We now show that £(■) is a fluid model solution on [0, L]. We have already 
verified property (i) of Definition 4.2. For any fixed n € N, 

sup (t { o},C(*)> < sup (l [0 , Kn) ,C(*)) 
te[o,L] te[o,L] 

< sup limsup(l [0iKn) ,C r (t)) 
te[a,L] f^oo 

(4.9) 

<hmsup sup (t[ , Kn ],C(t)) 
f-*oo te[o,L] 



where the second inequality follows by (4.3) and the Portmanteau theorem 
([1], Theorem 2.1). The last inequality follows by Lemma 3.10, since n(r) — > 
oo, as r — > oo, which implies that B r n ^ C B T n for all sufficiently large f. Since 
n G N is arbitrary and e n J. as n — ► oo, this proves property (ii) of Definition 
4.2. It remains to verify properties (iii) and (iv) there. We must show that 
either £(■) = on [0,L], or that for all t € [0,L], ((*) ^ and (4.2) holds. 
Since we have shown that C( - )) is constant on [0, L] and that property (ii) 
of Definition 4.2 holds, it suffices to show that (4.2) holds when £(0) ^ 0. 
For this, suppose that £(0) 7^ and note that property (ii) of Definition 
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4.2 implies that (x,((t)) = (x,C(0)) > for all t £ [0,L]. So since ((■) is 
continuous, 

(4.10) inf (l,C(t))>0. 

te[o,L] 

To show (4.2), we follow the analogous proof in ([6], Section 5.3) closely. We 
first restrict our attention to the case g £ V (see Definition 3.7) and derive 
a prelimit version of (4.2) which is satisfied by C r (') f° r an sufficiently large 
r. We will then pass to the limit in this relation to obtain (4.2) for £(•) and 
g £ V. Finally, we make a simple extension from V to C 

Let g £ V and note that since V C C, g' is uniformly continuous on M [recall 
that we extend g to be identically zero on (— oo, 0)]. So, there is a continuous 
nondecreasing function tp g : [0,oo) — ► [0, oo) with ip g (0) = 0, such that for 
any h € K, 

(4.11) sup||(/(a; + /i)-(/(x)||<^(H). 

Recall that for each f £ 7t, C r (') is a realization of fl r ' m (-) for some m < |/^J 
and some cu £ B^-y Since the sequence {C (')}fe7l remains fixed for the re- 
mainder of this proof, it is understood that all random objects indexed by f £ 
1Z and m < \ fT\ are evaluated at the particular oj £ B 1 ^,--, and m for which 

C(-) is a realization. Fix i £ [0, L] . For any N £ N and j £ {0, 1, . . . , N - 1}, 
define tj = and P =tj + \. Then we have for each f £ 1Z, 

N-1 

(g,C(t)) - (g,C(o)) = E ((9,C(t j )) - (g,?^))) 

3=0 
N-1 

= T,((9,C^))-(g(--S^),C(t J ))) 
j=0 

N-1 

(4-12) + J2((g(--S^),C(t j ))-{g,C(t J ))) 

j=0 

N-1 1 fE f ' m {P) 



E~ E 

3=0 



r UTf- 1 -m,P' 



N-1 
3=0 

where the first term in the last equality is by (3.2) (using t = tj, h = t 3 ' — tj 
there) and the fact that (1(0,00)5) = <?> since g(0) = for g £ C. By (4.3), 
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(4.13) lim S up||(l,C r ~(-)>llL<^T, 



(4.14) 



limsup||(l,C r X-))" 1 |lL<o°- 



By (4.13) and (4.14), we can assume for the remainder of the proof that r 
is large enough so that for some > 0, 



(4.15) 



sup (l,C f ( 5 ))V(l,C f (s))- 1 <M c . 
se[o,L] 



We handle the two terms in (4.12) separately. To begin with, since g E 
V C C has been extended to be an element of C^(R), we have the following 
first-order Taylor expansion for each j E {0, 1, . . . , N — 1} and each x E R+: 



(4.16) 



where hj = — S t ' tj and E R is in the interval [x — S t ] fj , x]. Note that by 
(2.20) and (4.15), 

(4.17) max \hA= max |^-™.| < Ih/I.^mx-Im < 

V ' j=0,...,N-l ' 31 j=0,...,N-l 1 *i-* j| -iV IIX ' S W/ " L - AT 

For each j E {0, ...,N — 1}, let Zj = sup se [ t C r ( s )) an d define hj = 



-Zjjj. Then 



N-l 



N-l 



J'=0 j=0 



h N % 



N-l 



(4.18) 



3=0 
N-l 



3=0 



5, 



n,/ 



For each iV E N and s E [0,t), let /jv(s) = Y^=q Zjt [tj p \{s) and define 
/iv(*) — 0- We can make the following estimate for the second term in (4.12): 

N-l N-l 

E - s£5) - <?(-u%)> - E c%)> 

3=0 3=0 



N-l 



< E sup \g(x - S r t 



g(x)-g'(x)h j \(l,C(t j )), 



j=0 xek + 
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N-l 

= J2 sup \g'{w*)h 3 -g'(x)h 3 \(l,C(t 3 )) 
N-l 

(4.19) <||(i,CT))ll L E wpMW-JWWthl 

+ W{x)\\hj - hj\), 

= M C (^(^)*M C 

+ iis'iico [f /«(») * - jTa.r w)- 1 *) V 

In the third line above we have used the Taylor expansion (4.16). The last 
inequality then follows from (4.15), (4.11), (4.17) and (4.18). The substi- 
tution of the second integral in the last line follows by (2.20) and (4.15). 
Let N — > oo in the above inequality. By the continuity of tp g and the fact 
that ip g (0) = 0, we see that the first term in the outer parentheses tends to 
zero. Note that /jv(s) — > (l,£ r (s)} -1 , as N —>■ oo, for any s £ [0, t) at which 
(1,(7* (-)) -1 1S continuous. Since it is continuous for almost every s [the path 
C( - ) is r.c.l.L], the second term in the outer parantheses tends to zero by 
(4.15) and bounded convergence. Furthermore, we note that 

N-l N-l , 

3=0 3=0 
and that as N — > oo, 

also by (4.15) and bounded convergence, since the function (g' , £"*(•)) (1, C r (0) 
is also continuous for almost every s. Together with the estimate (4.19), this 
implies that as TV — > oo, 

(4-20) J2(g(- - §1;™) - g(.), CM) -> - jT 

We handle the first term of (4.12) in a similar (although simpler) fashion. 
Once again we can use a first-order Taylor expansion for each summand 
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appearing in this term: 



(4.21) 
where h\ 



g{vl - s. 



mf- 



-m,P < 



- S '? __ i . , , and w\ E \v r - — SI 

U:r 1 —m,P 7 • 



• 3 - "uff-i-m,* ' ""J - l"« ~Ufr-i-m,P ' <1" SillCe ^ ~ ^ * ' 

w)| < i/-^ f° r each pair j, i in the first term of (4.12), we have by (2.20) and 
(4.15) as before that 

"3 1 



(4.22) 



maxl^lK^-IKl,^-))- 1 ! 



r < 
~ N 



Using the Taylor expansion (4.21) along with (4.22), we have 



'N-l 



rE f < m (P) 



0=0 i=f.E f . m (t,)-|-l 



5« - 3 



Uff' 1 -m,P' 



fE f ' m (t) 

E 

2.=r\E r > m (0)+l 



(4.23) 



JV-l 



fE f > m (P) 



E 5 E s'K)^ 



i=0 i=fE f . m (tj)+l 



C 



which tends to zero as N — > oo. By combining (4.20) and (4.23) above, we 
can let N — > oo in (4.12) to obtain the relation 



(4.24) (g,c(t)) = (9,em 



(t) 



o <u f oo> 



i=rE f - m (0)+l 

for all sufficiently large r. We would like to let r — ► oo in this relation to 
obtain 



(4.25) 



(g,C(t)) = (g,ao)) 



■ ds + at(g, v). 



lo (l,((s)) 

By (4.3), we have that the left-hand side, as well as the first term on the 
right-hand side of (4.24), converges to the corresponding term in (4.25). 
Similarly, (4.3), (4.13), (4.14) and (4.10) imply that the integrands in the 
second term on the right-hand side of (4.24) are uniformly bounded, and 
converge pointwise on [0, t] to the integrand in the second term on the right- 
hand side of (4.25). Thus the integrals converge by bounded convergence. 
To see that the third term on the right-hand side converges, note that since 
g E V, g = g)*. for some k E N (see Definition 3.7). So for any r ETZ large 
enough so that n(r) > k, 

fE f ' m {t) 

E 0fe( v i)-o*(0fe> 1/ ) 
i=fE f v' l (0)+l 
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< 



1 

r 

+ 



fE r - m ( L*/'n(f ) J l n (r) ) 

9k ( v i) -ult/ln(r)]ln(r)(9k,v) 
~ X) Sfc «)-«(*- L*An(f)J^n(f))(5fe,^> 

i=^' m (L*/'u(f)J'n(f)) + l 



(4.26) < E n[f) + ||^|U((^(t) - ^(L^n(f)J^n(r))) + <(r)) 

<en(f) + lbfcl|oof sup (E^ m (t + l n(f) )-E f ' m (t)) + al n{f) 



te[o,L] 

^ e n(f) + llSjfe l|oo£fi(f)) 

where the second inequality is by (3.75), and the last inequality follows 
[after relaxing the bound to £ n r?)] by (3.74) and (3.73). Since (4.26) holds 
for all f G 1Z sufficiently large, and since £ n rp) I as f — > oo, we obtain (4.25) 
from (4.24) by letting r — ► oo. 

We have shown that if C(0) ^ 0, then (4.2) holds for all t G [0, L] and g G V. 
We now extend to the case g G C, using the case g G C as an intermediate 
step. Suppose g G C and choose functions ^6V, k G N, such that {gk}kLi 
and {fi^}^! are uniformly bounded, and such that as k — > oo, 

(4.27) 5 fc — and g' k — >g', 

pointwise on (see Definition 3.7). We have established that for each k G N 
and t G [0,L], 

(4.28) (g fc , C(t)) = (9k, C(0)) ~ /' ( ,f '5/^ ^ + at(g fc ,^). 

It then follows by (4.27) and bounded convergence that, as k — > oo, the 
left-hand side as well as the first and third terms on the right-hand side of 
(4.28) converge respectively to the corresponding terms of (4.2). Similarly, 
since £(•) is continuous and £(t) ^ for all t G [0, L], the integral term also 
converges by (4.27) and two applications of bounded convergence. Thus, we 
obtain (4.2) from (4.28) by letting A; ^ oo. 

Finally, fix g G C. For k G N, choose a function tpk G C£(R+) such that 
ipk(x) ^ [0, 1] and — 2 for all x G K + , Vfc = 1 on [0, k] and tpk = on 

[A; + 1, oo). Let = tpkg, and note that {gk}kLi an d {<7fc}fcLi are uniformly 
bounded and that gk — ► g, and <^ — > g' are pointwise on R + as A; — > oo. 
Since has compact support, we see that gk G C for each fc G N. Therefore, 
(4.2) holds for all t G [0,L] and for each of the functions gk- By the same 
argument as that appearing after (4.28), this implies that (4.2) also holds 
for g. This completes the proof. □ 
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The next lemma asserts that for large r, the fluid scaled processes p> r ' m (-) 
can be uniformly approximated on [0, L] by paths in This fact, together 
with a result due to Puha and Williams [16] about the uniform convergence 
to steady state for elements of (see Proposition 4.5) constitutes the 
second main ingredient for proving state space collapse. 

Lemma 4.4. Assume (Q.l) and let T, L > 1 and < rj < 1 be given. Let 
{e n }^ =1 C (0, 1) be a sequence such that e n j 0, as n — > oo. Let l n , p, Mt, K n , 
7 n; N n , r n be the constants, and {B^} be the events given by Lemma 3.8. 
Then for each n € N, there exists r' n > r n such that for every r > r' n and 
C(-) £ W L , there exists £(■) 6 33l such that 

(4-29) ||d[C r (-),C(-)]llL<£n, 

(4-30) ||<x,C r (-)>-<x J C(-)>IU<^- 

Proof. We follow the proof of ([2], Lemma 4.1). Suppose on the con- 
trary that there exists n € N, a subsequence IZaTZ, and a sequence {( r (■)}- C 
Dl(-A/1f), such that ( r (-) G £$ r L for each f, and such that either 

inf inf ||d[C f (0,C(0]lli>^n or 

(4-31) 

inf inf |Kx,C r (0)-(x,C(-))ll i >^n. 

By Corollary 3.16 and Definition 4.1, there exists a £(•) € and a further 
subsequence {fj}j G ^ C 1Z such that as j — > oo, 

(4.32) C^'(-) ^CO- 

By Lemma 4.3, £(■) is continuous. So as j — > oo, 

(4-33) || d [^(.),C(-)]|| L -0, 

which implies by Lemma 3.5 and (3.76) that for each t € [0, L], 

(4.34) |( x , C ^))_( XjC (t))|^ asj^oc. 

Since (x>CG)) is constant by (ii) of Lemma 4.3, we have by (4.34) and (3.81) 
that in fact, 

(4.35) || <X, C^" (0> - <X, C(0> Hi — asj-oo. 
Together, (4.33) and (4.35) contradict (4.31). □ 

We now apply a result obtained in [16], on the rate at which fluid model 
solutions converge to their steady state, to choose L > 1 large enough so 
that elements of are nearly in steady state for t G [L — 1,L]. We first 



52 



H. C. GROMOLL 



recall the definition of the mapping S : Af F — ► D qo (A^f) introduced in [6], 
Lemma 4.9. Recall that the set Afp C A^f consists of all nonatomic £ G AIf, 
that is, all £ € Mf satisfying (1{ X },£) = 0, for all x G E + . Given £ G Alp, 
we define H(£)(-) = Cg(-), where Q(') S D 00 (A1f) is the unique fluid model 
solution for critical data (a, v) (defined on [0, oo) as in [6], Section 3.1), such 
that Cf (0) = £. The following proposition is a direct consequence of Theorem 
1.3 in [16]. Recall that the lifting map A y was defined in Definition 2.2. 

Proposition 4.5. Assume (Q.l) and let T > 1 and < 77 < 1 be given. 
Let {e n }^ =1 C (0, 1) be a sequence such that e n J, 0, as 00. TTien /or any 
e > 0, i/iere exists L* > 1 swc/i that for all £(•) G i G [L* — 1,L*] implies 

(4.36) d[C(t),A„<x,C(*))]<e. 

Proof. Let Z re , p, My, 7 n , iV n , r n be the constants, and {B r n } be 
the events given by Lemma 3.8. Let q and M q be the constants given by 
Lemma 4.3, and define 

(4.37) S = {£ G M C F : (1, £> V < x , £> V < MJ. 

Let ^ C D 00 (A / Jf) be the set of all fluid model solutions £(•) for critical 
data (a,v) (defined on [0, 00) as in [6], Section 3.1) such that C(0) G S. Then 
as a direct consequence of Theorem 1.3(i) in [16] (with M = M q and e = q 
there), we have that given 5 > 0, there exists L* > 1 such that 

(4.38) sup p(C(t),A„<x»C(0)» <6 for all t>L* — l, 

where /?(•,•) is the Prohorov metric on AJf defined in [16], Section 1. Note 
that for any £(•) G ^, A^(x,C(0)) G M«, where 

(4.39) M« = {?€M„:( X ,0<M ? }. 

Thus, to replace 5 and /?(•, •) with the given e and the metric d[-, •] in (4.38) 
above, it suffices to show that there exists 5 > such that for all £ G M£, 
the p-ball B p (£,5) of radius 5 centered at £, is contained in the d-ball 
^d(£> £ ) of radius e centered at £. To this end, let £ 9 = M q (x,v e )~ l v e be 
the element of with the greatest mass. Note that the metrics p(-, •) and 
d[-, •] both induce the same topology on AIf- So, we can choose 5 > such 
that B p (^ q ,5) C B d (£ q ,e). Suppose that £ G M£ and p((,£) < 5 for some 
C G Mf- Since M<? = {cv e : c G [0, M,(x, ^e)" 1 ]}, there is a £' G M$ such that 
£ + £' = £ 9 , where for two elements £1, £2 G A4f, we define £1 + £2 G AIf by 
(6 + 6) (-4) = + £2 (-4), for any Borel set Acl + . Using the definition 
of the metric ■), it is not difficult to verify that 

p(c+e,c+e)<p(c,o<s, 
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which implies that ( + £' G B p (£ q ,6) C B d (£ q ,e). By the definition (1.4) of 
d[-, •], this yields 

d[c,£] = d[C + £',£ + £']<£• 

Thus, we can choose L* > 1 so that 

(4.40) sup d[C(t),A„(x,C(0))]<e for all t>L* — l. 

C(-)G* 

Let be given by Definition 4.1 of this section. In [6], it was shown that 
for a fluid model solution £ (•) G D 00 (A^f) ) C(0 nas no atoms for all i > ([6], 
Lemma 4.3 and equation (4.33) ff.). This result carries over to fluid model 
solutions on [0, L*\ in a straightforward manner. So by Lemma 4.3 of this 
section, we have for any ((■) G that G S C Alp for all t G [0,L*]. 

For any C(-) G ^x*, let C(-) G D 00 (A^ F ) be defined by 

/CW, *e[o,L*], 

\H(C(L*))(t-L*), te(L*,oo). 

Note that if ^(0) = 0, then by Lemma 4.3, Definition 4.2 and [6], Theorem 
3.1, C(t) = for all t > 0. If C(0) ^ 0, then £(•) satisfies (4.2) for all t G [0, L*] 
and H(C(L*))(-) satisfies (4.2) for all t > 0. It is not difficult to see that in this 
case C(-) must also satisfy (4.2) for all t > 0. So, C(-) is a fluid model solution 
for critical data (a, z^) on [0, oo) such that £(t) = C(i) for all t G [0, L*]. Since 
C(0) G S for all CO) € @ L *, we have CO € & for all CO) G =^L*- Thus, by 
Lemma 4.3(h) and (4.40), for any £ G [L* — 1,L*], 

(4.42) sup d[C(t),A„<x,C(t)>]<e. 

C(.)P.«r. 



(4.4i) at) 



4.2. 5iaie space collapse. Before proceeding to the proof of state space 
collapse, we will need the following technical lemma, which provides a uni- 
form continuity property for the mapping H:Aip — ► D oc (A / Jf) on a set 
of measures £ G Alp that are close to the truncated invariant manifold M£. 
Recall that for any £ G Alp, = is the unique fluid model solution 

defined on [0, oo) such that C?(0) = £■ 

Lemma 4.6. Assume (Q.l) and ZeZ T,L>1 and < n < 1 &e given. Let 
{e n }^ =l C (0, 1) oe a sequence such that e n [ 0, as n — > oo. Lei Z n , p, Mr, K n , 
7nj r n ^ e ZTie constants, and {B^} be the events given by Lemma 3.8, 
and let q and M q > 1 be the constants given by Lemma 4.3. Let M\ > M q 
and define 

(4.43) M* = {eeM„:<x,0<Mi}. 

TZien /or any e > 0, Z/iere exists (5 > sucZi that for all CO) G and £ G 
satisfying 

(4.44) d[((0U] v | ( x ,C(o))-(x, 
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we have 

(4.45) sup d[((t),Q(t)]<e. 

te[o,L] 

Proof. Fix e > and choose no such that e no < e/4. Let 

(4.46) §o = s(K no M)(( X ,u e )Al) _ 

8 

Note that by (4.9), any ((■) G ^ L satisfies 

(4-47) sup (1 ( >,£(*))<£=&<£. 

te[o,L] 2 £5 

Note also that for any ((■) G C(<) = s (C(0))(i) for all t G [0,L]. By The- 
orem 3.8 in [15], S is continuous on .M F . Let £i = Mi(x, v^)~ x v e . Since 
£l G .Mp, we can choose < 5\ < 1 such that for any £ G ,M F satisfying 
d[£,£i] < we have 

(4.48) sup d[Ce(t), (*)]<£■ 

tG[0,2MiL/5 ] 

Define 

(4.49) 5 = Ml, 
v ; 2Mi ' 

and consider any £(■) G such that for some £ G M*, (4.44) holds. The 
argument splits into two cases. First, suppose that (x>C(0)) < <5o- Then, by 
(4.44), (4.49) and the fact that Mi > 1, (x,£) < 2<5 . This implies by (4.46) 
that (1,0 = ( X ,0(x,^y 1 <e/A. Thus, 

sup d[C(*),Cf(t)] 
te[o,L] 

< sup d[C(*),0]+d[0,e] 

te[o,L] 

(4.50) < sup 2(1, C(*))+ 2(1,0 

te[o,L] 

<2 sup ((l[o,«„ ),C(*)} + (l(« no> oc),C(*))) + ^ 
te[o,L] z 

<| + — (x,CW) + £ 



4 k„ x 2 



<£, 



where the first inequality uses the fact that Ct(t) = Q(0) = £ f° r an £ > 0, 
since £ is on the invariant manifold ([16], Theorem 1.1). The second inequal- 
ity uses definition (1.4), the fourth uses (4.47) and Markov's inequality, and 
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the last inequality uses (4.46) and the fact that since £(■) is a fluid model 
solution, (x,C(t)) = (XiC(O)) < for all t 6 [0,L] ([6], Theorem 3.1). This 
establishes (4.45) in the first case. 

For the second case, suppose that (x, C(0)) > So- Since 6\ < 1 < Mi, (4.44) 
and (4.49) imply that (x,Q > <5 /2. So £^ 0. Let c x = Mi(x,£) _1 and denne 
the transformed function £(•) : [0, ciL] — > Mf by 

(4.51) C(t) = ciC(*/ci), t€[0, Cl L]. 

Using a change of variables, it is not difficult to verify that £(•) satisfies (4.2) 
for all t G [0, ciL] and is therefore a fluid model solution on [0,ciL]. Note 
that C(0) = ciC(O), and that ci£ = £i. Also note that since £ G M*, we have 
ci > 1. Thus, using definition (1.4), we have 

d[C(0),£i] = d[ Cl C(0), Cl £] 

, . < Cl d[C(0),£] 
4.52 

A 

do 
= St. 

This implies that 

sup d[c(t),c*(t)]< d[ciC(o,cie] 

t£[0,L] te[o,L] 

(4.53) = sup d[C(t),Cfr(*)] 

te[o,ciL] 

<£■ 

The first line in the above inequality uses (1.4), the fact that c\ > 1 and the 
fact that Q(i) = £, since £ € M*. The second line follows by definition of 
CO), and the fact that ci£ = £i G M^. The last line then follows by (4.52) 
and (4.48) [with £ = C(0) there] by noting that c x L < 2M 1 L/6 . □ 

We are now ready to prove state space collapse for the sequence of dif- 
fusion scaled state descriptors {/i r (-)}- This leads directly to the proof of 
Theorem 2.3, which appears at the end of the section. 

Theorem 4.7 (State space collapse). Assume (Q.l) and let T > 1 be 
given. Then as r — > oo, 

(4.54) ||d[/rQ,A^(-)]|| T ^0. 

Proof. We must show that for any < £*,rj < 1, there exists r* E 1Z 
such that r > r* implies 

(4.55) P'-(||d^(-) > A^(.)]|| T <0>l-^ 
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Fix < £*, r] < 1. We first note that by definition of the metric d[-, •] and the 
lifting map A„, 

(4.56) d[A u wi,A u w 2 ] < c u \wi - u>2\, w\,W2 G M+, 

for some constant c u > 1. Choose a sequence {e n }J^ =1 such that e n j as 
ra — > oo. By Proposition 4.5, using e = e*/3 there, we can choose L* > 1 such 
that for all C(-) G ^i*, i G [L* - 1,L*] implies 

(4.57) d[((t),A u (x,((t))]<j- 

Let Z n , p, My, K n , j n , N n , r n be the constants, and {B^} be the events given 
by Lemma 3.8. Let q and M q > 1 be the constants given by Lemma 4.3, and 
let 5 > be given by Lemma 4.6 for L = L* , Mi = 2M g and e = e*/l2. Fix 
n* G N large enough so that 

(4.58) e n . < - A . 

Choose r* large enough so that n(r*) > n* and r* > r' n ,, where r' n , > r n * is 
given by Lemma 4.4. Fix r > r* for the remainder of the proof. The argument 
splits into two cases. We first consider the process /x r '°(-) on B^\, and then 
consider the processes /2 r,m (-) for m € {1, ... , L r ^J} on ^n(rV 

For the first case, consider any realization C r ( - ) G of the process p> r '°(-) 
on the event B r n ^ r y By (3.77), there exists a £ € Mj, such that 

(4.59) d[C r (0),£] V |(x,C r (0)) - (X,0l <e»(r) <£»*• 

By (3.76) and the fact that e n _* < 1 < M T < M q < M u we see that £ G = 
{£ G Mj, : (x,0 < 2M q }. Let Q(-) = =(£)(■) and note that since £ is on the 
invariant manifold, £ = ^(0) = Q(t) f° r all ^ £ [0, £*] ([16], Theorem 1.1). 
So by definition of M. u and A v , 

(4.60) Q(i) = A„(x,£(t)) for all t G [0,L*]. 

Let £(•) G =^L* be given by Lemma 4.4 for C(-) an d n = n*. Note that since 
Q(-) and £(•) are fluid model solutions on [0,L*], we have for all t G [0, L*] 
([6], Theorem 3.1), 

(4.61) (x,Cc(t)) = (x,^(0)> and (x,C(*)> = <X,C(0)>. 

By combining the result of Lemma 4.4 cited above with (4.59) and (4.58), 
one obtains 

d[c(o),e]<d[C(o),c r (o)]+d[C(o),£] 

(4.62) < 2e n * 

<S, 
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l<x,Q(0)>-(x,C(0)>| 

tA , < \(x,cm) - (x,c r (o))\ + \(x,cm - (x,c(o))i 

(4.63) 

<2e n . 
<S. 

This implies, by choice of <5 and Lemma 4.6, that for any t 6 [0, L*], 
(4-64) d[m ,Q(t)}<^. 

Since r > r*, we have for any t G [0, L*], 
d[C(t),A v { X ,C(t))] 

< d[C(t), at)] + d[C(t), Q(t)\ + d[Q W, A,(x, Cc(*)>] 

, . + d[A v <x, &(*)>, A„<x, C(*))] + d[A,( x , C(i)), A„( X , C r W>] 

e* 

- £n * + 12 + ° + c ^ 2£n * + Ci/£n * 
e* 

where the first estimate in the second inequality above follows by Lemma 
4.4, the second estimate follows by (4.64), the third by (4.60), the fourth by 
(4.56), (4.61) and (4.63), and the fifth estimate uses (4.56) and Lemma 4.4. 
The last inequality follows by (4.58) and the fact that c„ > 1. 

We now proceed to the second case. Let C r (') € SB r L » be a realization of 
the process /U r,m (-) on B^,y where m G {1, [r^j}- Let £(•) € be 
given by Lemma 4.4 for C r ( - ) an d n = n* . Since r > r*, we have for any 
te[L*-l,L*}, 

d[c r (t),A u ( x ,C(t))] 

<d[c r (t),C(t)]+d[C(t),A,( x ,c(t))] 

(4.66) +d[A u ( x ,at)),A u (xX r (t))] 

e* 

< £n* + — + C v E n * 

<£!, 

~ 2 

where the three estimates in the second inequality are by Lemma 4.4, (4.57), 
(4.56) and Lemma 4.4, and the last inequality is by (4.58) and the fact that 
C>1. 
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Finally, combining the estimates (4.65) and (4.66), we have on BWs, for 
r > r*, 

||d[/T(-),A^(.)]|| r 

= sup d[fl r (t),A u (x,fl r (t))] 
te[0,rT] 

(4.67) < sup d[^°(t),A u ( X ,fl^(t))} 

te[o,L*] 

+ sup d[ff> m (t),A u ( X ,^ m (t))} 

l<m<[rTJ, t€[L*-l,L*] 

e* e* 

< — I — . 

- 2 2 

Since P r (B 7 n ^) > 1 -rj for r > r* , (4.55) is proved. □ 



Proof of Theorem 2.3. By Proposition 3.1, we have that W r {-) => 
W*(-) as r — > oo, where W*(-) is a reflected Brownian motion in ~R + with 
drift —A, variance aa 2 + /36 2 and initial condition VF*(0) equal in distribu- 
tion to {x, O) ■ Since A v : M + — ► A^f is continuous, the continuous mapping 
theorem implies that A u W r (-) =^ = A U W*(-) as r — > oo. Thus, The- 
orem 4.7, combined with the "converging together lemma" ([1], Theorem 
4.1), implies that jl r {-) as r — > oo. □ 
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